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Introduction 


The deep space environment is a very hostile one for lifd as it exists on 
earth. The protection of space travelers from this hostile environment is a 
challenge to biology and technology in which both disciplines are dealing in 
largely uncharted territory. Solutions require communication between these 
disciplines of unprecedented effectiveness. In technology, the analysis of 
the behaviour of complex systems is a well-developed technique which yields 
many benefits. An important benefit is the ability to simulate the behaviour 
of a complex system under almost any imaginable set of circumstances using 
mathematical models with great savings in time and expense. The accuracy 
of such models of technological systems can be very high since the character- 
istics of their components are well known, and the controlling systems have 
known structure and characteristics. 

It is highly desirable to develop similar mathematical models of biological 
systems for obvious reasons: simulations of the biological responses to 
dynamically varying sets of environmental circumstances would allow savings 
in time and expense required for evaluating the effects of such environments. 
Another very attractive feature of mathematical models of biological systems is 
that they can be made to communicate with models of technical systems, and 
it thus becomes possible to simulate the combined system of man with his 
life support system. 

There are serious difficulties in constructing mathematical models of biological 
systems. These difficulties rest in the complexity of biological systems with 


respect to the number of identifiable components, but they are especially 
due to a great deal of redundancy in the control systems. This redundancy 
results in control systems with many multiple control loops. In a multiple 
control loop system, it is hard to study the system by opening one or more 
control loops since one cannot be sure that there are no remaining loops and 
that a particular component under study is really isolated. 

It is thus not surprising that a mathematical model of thermoregulation of 
man, such as is described in this report, is far from accounting completely 
for all thermoregulatory responses ever reported. It can, however, give 
reasonable estimates of dynamic thermophysiological responses to a variety 
of environmental conditions and rates of internal heat production. This 
report is intended to present the current status together with the experimental 
findings which form the basis of much of the structure of the current model. 

Development of the Model 

An early forerunner of the current model was described in 1966 (1). This 
forerunner was implemented on an analog computer, but as the definition of 
the model increased, it soon became preferable to express the model equations 
in Fortran so that the simulations could be carried out on a digital computer. 

As proposed for models of the respiratory system by Yamamoto and Raub (2), 
symbol tables and equations will be given in a Fortran notation. 

Although in practice the two are very closely interwoven, we can conceptually 
distinguish two separate systems in thermoregulation: the control system and 
the controlled system. The controlled system is a representation of the body 


2 



in the thermal characteristics of its different parts. Thermal stresses act on 
the controlled system and cause strains. These thermal changes are recognized 
by sensor mechanisms belonging to the control system and after integration 
the control system causes corrective actions which act to reduce the thermal 
strain. 

The controlled system in the current model consists of six segments, each 
with four concentric layers. A central blood compartment exchanges heat with 
all other compartments via the convective heat transfer occurring with the blood 
flow to each compartment. Each of the twenty-five compartments is represented 
by a heat balance equation which accounts for conductive heat exchange with 
adjacent compartments, metabolic heat production, convective heat exchange 
with the central blood compartment and evaporative heat loss and heat exchange 
with the environment if the compartment is in direct contact with the environment. 

The Controlled System 

The controlled system is based on a man with a body weight of 74.1 kg and 
2 

a surface area of 1.89 m . It consists of the head which is considered as a 
sphere, and cylinders representing the trunk, the arms, the hands, the legs and 
the feet. Both arms, hands, legs and feet are represented by one cylinder each. 
Table 1 gives a list of symbols used in the controlled system with their definition 
and dimension. 
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TABLE 1 


LIST OF SYMBOLS USED IN THE CONTROLLED SYSTEM WITH DEFINITION AND 


DIMENSIONS 


Sym- Vector 
Ppl_ Length 

Definition 

Dimensions 

c (N) 

25 

Heat capacitance of compartment N 

kcal . °C _1 

T (N) 

25 

Temperature of N 

°C 

°C . h _1 

F (N) 

25 

Rate of change of temperature in N 

HF (N) 

25 

Rate of heat flow into or from N 

kcal . h 1 

, -1 -1 

TC (N) 

24 

Thermal conductance between N and N + 1 

kcal . h . °C 

TD (N) 

24 

Conductive heat transfer between N and N+l 

kcal . h" 1 

QB (N) 

24 

Basal metabolic heat production in N 

kcal . h 1 
kcal . h 1 

Q (N) 

24 

Total metabolic heat production in N 

EB (N) 

24 

Basal evaporative heat loss from N 

kcal . h _1 
kcal . h * 

E (N) 

24 

Total evaporative heat loss from N 

BFB (N) 

24 

Basal effective blood flow to N 

1 . h _1 
1 . h" 1 

BF (N) 

24 

Total effective blood flow to N 

BC (N) 

24 

Convective heat transfer between central 
blood and N 

kcal . h" 1 

HC (I) 

6 

Convective and conductive heat transfer 
coefficient for Segment I 

kcal . m ^ . h \ °C 1 

S (I) 

6 

Surface area of Segment I 

2 

m 

, -2 -1 -1 

HR (I) 

6 

Radiant heat transfer coeff . for Segment I 

kcal . m . h . °C 

H (I)' 

6 

Total environmental heat transfer coeff. 
for Segment I 

kcal . h _1 . °C _1 
-1 

m . sec 

V 


Air velocity 

TAIR 

RH 


Effective environmental temperature 
Relative humidity in environment 

°C 

TIME 


Elapsed time 

h 

PAIR 


Vapor pressure in environment 

mm Hg 

ITIME 


Elapsed time 

min 

INT 


Interval between outputs 

min 

DT 


Integration step 

h 

P (I) 

10 

Vapor pressure table from 5 - 50°C 

mm Hg 

EMAX (I) 

6 

Calc. max. rate of evaporative heat loss 
from Segment I 

kcal . h 1 
kcal . h _1 

WORK 


Total metabolic rate required by exercise 

pskin(i) 

6 

Saturated water vapor pressure at skin temp. 

mm Hg 
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A schematic representation of one segment is given in Figure 1, showing the 

interrelations between the four concentric layers. The subscript I in Table 1 

and' hereafter refers to the segment, 1=1 being the head, 1 = 2 the trunk, 1 = 3 

the arms, 1 = 4 the hands, 1 = 5 the legs and 1 = 6 the feet. The subscript N is 

* 

used to indicate the individual compartments, such that N = 4 1-3 always 

* 

identifies the core layer of segment I, with N = 4 1-2 representing the muscle 

* * 
layer, N = 4 I - 1 the subcutaneous fat layer and N = 4 I is the skin com- 
partment of segment I. The central blood compartment is represented by N = 25. 



Figure 1. Schematic diagram of the four compartments of Segment I. 
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The first step in the development of the characteristics of the various 
segments is the determination of their relative surface areas and volume. 
Surface areas of men and women have been measured by many investigators, 
and Table 2 gives such values as selected from Meeh (3) and DuBois and 
DuBois (4) for men, and by Bradfield (5) for women as averages for 10 indi- 
viduals each. 


Table 2 

Values for Surface Areas of Men and Women 


Segment 


S u r f a 

o 

CD 

> 

CD 

a 

1 0 

men 

1 0 

women 

2 

m 

% of total 

2 

m 

% of total 

Head 

0.1326 

7.00 

0.1129 

6.49 

Trunk 

0.6804 

36.02 

0.6279 

36.07 

Arms 

0.2536 

13.41 

0.2210 

12.70 

Hands 

0.0946 

5.00 | 

0.0783 

4.50 

Legs 

0.5966 

31.74 

0.5904 

33.92 

Feet 

0.1299 

6.86 

0.1100 

6.32 

TOTAL 

1.8877 

100.00 

1.7404 

100.00 


Volume measurements of the various segments were not found in the literature in 
usable form, so that such measurements were approximated by measuring on five 
men and five women the circumference of the head, the length and two circumferences 
of the trunk, the length and circumference of the upper arm, and the length and two 
circumferences of the forearm. Similarly, the length of the thigh from crotch to knee, 
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and two circumferences were measured for the thigh, and a similar set of 
measurements were taken for the lower legs; the volume of hands and feet were 
measured by displacement upon submersion in water. The volumes of the other 
segments were estimated from lengths and diameter assuming cylindrical shapes. 
The total estimated volumes were checked against the individual weights and 
were found to agree to within 8% with the total weight assuming a density of 1.00. 
Table 3 gives the percentages of the total volume in the various segments, as 
averaged for 5 men and 5 women. 


Table 3 

Percent of Total Body Volume in Six Segments 


Segment 

Men 

Women 

Average 

Head 

5.34 

5.58 

5.46 

Trunk 

56.70 

57.00 

56.85 

Arms 

7.78 

6.52 

7.15 

Hands 

0.88 

0.90 

0.89 

Legs 

27.90 

28.70 

28.30 

Feet 

1.43 

1.34 

1.38 

TOTAL 

100.00 

100.00 

100.00 


According to Wilmer (6), the tissue distribution in the human body is 25% 
skin and fat, 11% viscera, nervous tissue 3%, muscle 43% and skeleton 18%. 

Scammon (7) tabulates average weights of different body divisions and organs 
and his numbers can be used to make assignments of the tissue types to the layers 
in each segment as given in Table 4. 
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TOTAL 5.00 15.00 43.00 15.00 8.00 14.0 37.00 100.00 












With the aid of the preceding tables, the necessary computations to obtain 
C (N) and TC (N) can be carried out for any size man; but, for our purposes, 
we will limit ourselves to a "Standard" man with a height of 172 cm and a weight 

2 

of 74.4 kg who, according to DuBois and DuBois (8) has a surface area of 1.89 m , 
The process for obtaining C (N) is indicated in Table 5. The values are based on 
a specific heat for the skeleton of 0,5 gcal/g . °C, 0.6 gcal/g . °C for fat and 
0,9 gcal/g . °C for all other tissues. The central blood compartment representing 
the blood in the heart and the great vessels is assumed to contain 2.5 liters of 
blood at a specific heat of 0.9 gcal/g . °C. This heat capacitance of 2.25 
Kcal/°C should be subtracted from the heat capacitance of the trunk core which 
thus becomes 9.82 Kcal/°C. 

Thermal conductance between layers in each segment can be calculated from 
the dimensions derived before, using the method outlined by Stolwijk and Hardy (1). 
The lengths of the segments were obtained from measurements taken for the volume 
estimates described above. Results are given in Table 6. Although the volume 
estimates and surface area estimates have different origins, there is a con- 
siderable consistency in the data as can be demonstrated by recalculating the 
total surface area from the estimates in length and radius given in Table 6. 

Ignoring end surfaces, such a calculation yields a value very close to the DuBois 
area assumed at the outset. 


9 


<0 .. 
o O 
O w o 



o S o 


O (0 

3 ° 

o 


U|o 

o 


°l u 


c So 

O O & o 



H 

CO 


00 

eg 

eg 

CM 


»— i 

o 

eg 

• 

• 

• 

• 

• 

• 

° 

r-H 

o 

o 

i-H 

o 


LO 

00 

CD 

o 


eg 

CO 


i-H 

eg 

eg 

• 

• 

• 

• 

• 

• 

o 

»— I 

o 

o 

i-H 

o 

eg 

LO 

00 

<D 

CO 

CO 

eg 

eg 

LO 

O 


i-H 

• 

• 

• 

• 

• 

• 

o 


o 

o 

rH 

o 

tv 



to 

00 

eg 

CO 

o 

<D 

r-H 

CO 

eg 

• 

• 

• 

• 

• 

• 

o 


O 

o 

eg 

o 

CO 

LO 


CD 

e*. 

CD ] 

CO 

i-H 

o 

o 

r-H 

O 

• 

• 

• 

• 

• 

• ' 

o 

CD 

CO 

o 

CD 

o 


I-H 





r^. 

O 

tv 

e* 

CD 

IV 

CO 

CT> 

CO 

o 

»-H 

o 

• 

• 

• 

• 

• 

• 

o 

t'v 

CO 

o 

o 

o 


1— 1 



i-H 


eg 

f-s. 

i-H 



CO 

eg 

o 


r-H 

eg 

eg 

• 


• 

• 

• 

• 

eg 

eg 

i-H 

o 


o 


r-H 






00 

CO 

CD 


CO 

o 

CD 

eg 

eg 

CD 


• 

• 

• 

• 

• 

• 

CO 


eg 

o 

to 

o 


r— i 





I-H 

CD 

CO 

CO 

CO 

LO 

to 

CD 

CD 

o 

IV 

o 

• 

• 

• 

• 

• 

• 

rH 

O 

o 

o 

i-H 

o 


i—t 





<D 

LO 


CO 

eg 

CD 


00 


o 

O') 

O 

• 

• 

• 

• 

• 

• 

I— 1 

rH 

o 

o 

i-H 

o 


i— 1 ( 





1— 1 

i— 1 

CO 

i-H 

rH 

00 

CD 


v- 

rH 

LO 

i-H 

• 

• 

• 

• 

• 

• 

O 

I-H 

o 

o 

eg 

o 

eg 

CO 

•—* 

CO 

eg 


eg 

00 

LO 

eg 

O 

CO 

• 

• 

• 

• 

• 

• 

rH 

eg 

r— 1 

o 

LO 

o 

eg 

o 

CD 


00 


o 

o 

O 

CD 

CD 

CD 

• 

• 

• 

• 

• 

• 


i-H 


O 

o 

o 





eg 



i — i 

LO 

CD 

00 

CO 

• 

• 

• 

• 

• 

• 

to 

CO 

CD 

o 

Cv 

t-H 


LO 



eg 





CO 



TJ 
f 0 
0) 

C 

s 

m 

TJ 
C 
f 0 

w 

S’ 

■M 

<D 

CD 

E 

H 

< 

E 


PH 






10 


Total 100.0 74.40 11.18 5.57 16.39 14.74 27.57 


































Table 6 

Calculation of Estimated Thermal Conductance Between Compartments 



Volume 

3 

cm 

Length 

cm 

— 

Radius 

cm 

Avg. cond. 
gcal/cm . 
°C . sec 
x 10 4 


Head 

core 

■a 


8.98 

10 

1.38 


muscle 

1 Imm 

— 

9.32 

10 

11.4 


fat 

B9 

— 

9.65 

8 

13.8 


skin 

4020 



9.88 

8 

— 

Trunk 

core 

14680 

mm 

9 

10 

1.37 


muscle 

32580 

■■ 

1 

10 

4.75 


fat 

39650 

60 

14.40 

8 

19.80 


skin 

41000 

60 

14.70 

8 

— 

Arms 

core 

mm 

m 

2.83 

10 

mm 


muscle 

■Mill 

■9 

4.48 

10 

9 


fat 

6580 

112 

4.85 

8 



skin 

7060 

1 

112 

5.02 

8 

mm 

Hands 

core 

mm 

96 

isa 

10 

5.50 


muscle 

m 

96 

mm 

10 

9.65 


fat 

480 

96 

1.27 

8 

9.90 


skin 

670 

96 

1.49 

8 

— 

Legs 

core 

6910 

160 

3.71 

10 

9.0 


muscle 

17100 

160 

5.85 

10 

12.4 


fat 

19480 

160 

6.23 

8 

64.0 


skin 

20680 

160 

6.42 

8 

— 

Feet 

core 

HI 

mm 

m 

10 

14.0 


muscle 


■B9 

Da 

10 

17.7 


fat • 

730 

125 

Da 

8 

14.1 


skin 

970 

125 

H 

8 
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The various segments have different dimensional characteristics and, as 
a result, the overall environmental heat transfer coefficient H (I) has different 
values. Radiant heat transfer depends on effective radiant surface areas which 
are high for head and trunk, somewhat lower for arms and legs and even lower 
for hands and feet. Estimates for the effective radiant heat transfer coefficients 
of the six segments are given in Table 7. It should be pointed out that these are 
estimates, and that behavioral adjustments can increase and decrease these co- 
efficients substantially. The values for the convective heat transfer coefficient 
given in Table 7 are based on measurements by Nishi and Gagge (9) in our 
laboratory using the naphthalene mass transfer method. They represent values 
in still air and we can assume that the effective air speed was 0.1m, sec * 
due to gravitational convection. The combined heat transfer coefficient H (I) 
can be calculated for any air speed V by the expression 

H (I) = (HR (I) + HC (I)* (V/0.1)** 0.5)* S (I) 

This expression is valid within about + 5% for an ambient temperature range 
of 10° to 40° C and for air velocities up to 10 m . sec 

The value for the combined heat transfer coefficient is slightly higher 
than found experimentally for sitting, resting subjects in shorts (10, 11), but 
this could be due to the minimal clothing and to behavioral adjustments of 
their effective surface area by the subjects. 
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The model will be called upon to compute heat loss caused by evaporation 
of secreted sweat. The amount of sweat secreted is determined by the con- 
trolling system. For any given set of conditions there is a maximum rate of 
evaporation which is possible, indicated by EMAX (I) for segment (I). This 
maximum rate depends on the vapor pressure in air PAIR, the vapor pressure 

* 

at the skin PSKIN (I) and the convective heat transfer coefficient H (I) - HR (I) S(I). 

The following expression will compute EMAX (I) for segment I: 

EMAX (I) = (PSKIN (I) - PAIR)* 2.2 * (H (I) - HR (I) * S (I) ) 

in which 2.2 is the Lewis relationship as described and validated by Brebner, et 
* 

al. (12). If E (4 I) for any segment I exceeds EMAX (I), then it must be reduced 

* 

to EMAX (I). The ratio E (4 I) /EMAX (I) for each segment gives the fractional 
wetted area, a factor first described by Gagge (13) and later implicated as an 
important determinant in thermal comfort by Gagge, et al. (14). A weighted 
average of this ratio will correspond to the wetted area described in these 
reports. 

Aschoff and Wever (15) have estimated that the metabolic rate of the 
brain is 16% of the total basal metabolic rate, and that the trunk core accounts 
for 56% of the resting metabolic rate. They estimate the total skin and mus- 
culature to produce 18% of the basal metabolic heat which leaves 10% for the 

skeleton and connective tissue. If the resting metabolic rate is presumed to 

-1 -2 -1 
be 74.4 Kcal . h for the whole standard man (39.4 Kcal . m . h ), the 

head core has a metabolic rate of 11.9 Kcal . h and the trunk core of 

41.6 Kcal . h *. A total resting metabolic rate of 13.4 Kcal . h 1 is to be 
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divided between skin and musculature. If the resting metabolic rate is 0.3 
Kcal . h 1 . kg 1 for skin tissue and subcutaneous fat is included in this, 
we find a total of 4.47 Kcal . h 1 to be divided between the various skin and 
fat layers. This leaves 8.93 Kcal . h 1 to be divided proportionally between 
the different muscle compartments. All core compartments can then share the 
remainder of the resting heat production in proportion to their weight. The 
values of QB (N) resulting from these estimates and considerations are given 
in Table 8. 

In the resting state, the total evaporative heat loss is about 18 Kcal . h 1 
(Stolwijk and Hardy (10) ). Of this, about half, or 9 Kcal . h 1 represents 
respiratory heat loss which is considered to come from the trunk core. The 
remaining 9 Kcal . h 1 comes from the skin and represents diffusion losses 
through the skin. Strictly speaking both of these values are affected to some 
extent by the water vapor pressure in the ambient air, but this effect is 
relatively small and is ignored for our present purposes. 

The values of basal evaporative heat loss EB (N) resulting from these 
estimates are given in Table 8. 

The controlled system is further characterized by convective heat exchange 
between its parts as a result of blood flow. Certain regions are characterized 
by a very constant blood flow. The brain receives a very constant flow of 
about 750 ml . mm 1 or 45 1 . h MiG). Resting muscle can be assumed to 
require 1.2 1 . h 1 of blood flow for each Kcal . h 1 of metabolic heat pro- 
duction, purely on the basis of supplying the needed oxygen. Stolwijk and 
Hardy (1) added up the blood flows to visceral organs in the trunk core and 
estimated a total flow of 210 1 . h *. The skin, in the resting state, receives 
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far more blood flow than is required to sustain its metabolic rate. Behnke 

and Willmon (17) found a total skin blood flow in resting man of 
-1 -2 

240 ml . mm . m at an ambient temperature of 3 5°C which would cause 

considerable vasodilatation. In the neutral environmental temperature of 

-1 -2 

30°C, this is assumed to be reduced to 105 ml . min . m in proportion 
to the reduction in core to skin conductance which occurs at an environ- 
mental temperature of 30°C in comparison with 3 5°C. Hertzman, et al, (18) 
investigated the distribution of blood flow to different skin areas, and skin 
blood flow was distributed in Table 8 according to his findings. 
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Table 8 

Estimate of Core Compartment Resting State Heat Production 
and Basal Evaporative Heat Losses 


Segment 

a) 

i 

Layer 

(N) 

Weight 

kg 

C (N) 

QB (N) 

1 

EB (N) 

1 

BFB (N) 

i 

Head 

core 

1 

1 

2.22 

12.84 


45.00 

(1) 

muscle 

2 

■ 

0.33 


— 

0. 12 


fat 

3 


0.22 

0.11 

— 

0.13 


skin 

4 

0.27 

0.24 

0.08 

0.63 

1.44 

Trunk 

core 

5 

■99 

9.82 

45.38 

9.00 

210.00 

(2) 

muscle 

6 

mum 

16.15 

5.00 

— 

6.00 


fat 

7 

BBS 

4.25 

2.13 

— 

2.56 


skin 

8 

■a 

1.21 

0.40 

3.25 

2.10 

Arms 

core 

H 

2.25 

1.41 

S9 

___ 

0.84 

(3) 

muscle 

mm 

3.37 

3.04 

■is 

— 

1.14 


fat 

li 

0.97 

0.58 

EES 

— 

0.20 


skin 

12 

0.48 

0.43 

0.13 

1.20 

0.50 

Hands 

core 

H 

mm 

m 

■ffi 


0.10 

(4) 

muscle 

■S 

■89 


■i 


0.24 


fat 

15 

0.15 

0.09 

189 


0.04 


skin 

16 

0. 19 

0. 17 

0.05 

0.45 

2.00 

Legs 

core 

17 

6.94 

4.24 

2.23 


2.69 

(5) 

muscle 

18 

10. 19 

9.17 

2.86 

— 

3.43 


fat 

19 

2.38 

1.43 

0.43 

— 

0.52 


skin 

20 

1.20 

1.08 

0.32 

2.85 

2.85 

Feet 

core 

21 

■H 

mm 

SI 

^ ^ 

0.16 

(6) 

muscle 

22 

wSm 

■US 

IMS 

— 

0.02 


fat 

23 

0.22 

0. 13 

0.04 

— 

0.05 


skin 

24 

0.24 

0.22 

0.07 

0.62 

3.00 


Central 








Blood 

25 

2.50 

2.25 

— 

— 

— 

Total 



74.45 

I 

59.17 

74.45 

18.00 

285.13 
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The Controlling System 


The controlling system in body temperature regulation can be divided 
into three distinctly separate parts. The first part contains the sensing mechanisms 
which recognize the thermal state of the controlled system. The second part 
receives information regarding the thermal state, integrates it and sends out 
appropriate effector commands to the various effector systems. The third part 
of the controlling system receives the effector commands and if appropriate 
modifies them according to the conditions at the periphery before translating 
such commands into effector action. 

Thermoreceptors have been described in a wide variety of locations in 
the body. Sometimes these descriptions are based on functional thermo sensitivity 
in short, local temperature changes and produce thermoregulatory effector output 
much in excess of the size of the stimulus. On other occasions, and using 
electrophysiological methods, investigators find fibers in which the electrical 
activity is determined by the temperature of tissues innervated by such a fiber. 

Areas in which thermosensitivity has been demonstrated by both methods 
are the preoptic area and anterior hypothalamus in the brain (19, 20), the spinal 
cord (21, 22) and the skin (23, 24), and the abdominal viscera (25, 26). Other 
areas have been suspected, such as muscle tissue or the veins draining muscle 
(27). Final and irrefutable proof of the functional role of thermosensitive areas, 
and the precise mode of detection is extremely difficult to obtain, even for the 
anterior hypothalamic area where there is general acceptance of the importance 
in thermoregulation. 
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Such uncertainties stem from the fact that it is usually difficult to dis- 
sociate temperatures in various parts of the body because of the strong thermal 
coupling which takes place. Uncoupling for experimental purposes is difficult 
because it can only be achieved for a very short time during rapid transients 
in load. During such transients, thermal fluxes are high, thermal gradients 
are steep, and the available time is short, all adding up to potential errors 
in such measurements and to general difficulties in devising appropriate 
methods. 

A listing of all symbols used in the description of the controlling system 
is given in Tables 9 and 10. 

In order to allow for maximum flexibility in the evaluation of different thermo- 
regulatory controller concepts, it is easiest to assume that thermoreceptor 
structures are present in all tissues. Error signals are derived from all compart- 
ments with the following expression: 

ERROR (N) = T (N) - TSET (N) + RATE (N)* F (N) 

In first approximation, the error signal is thus equal to the difference between 
the instantaneous temperature T (N) and the reference temperature TSET (N). 

For appropriate compartments where dynamic sensitivity of the thermoreceptors 

•k 

is quantitatively known, the dynamic term RATE (N) F (N) can assume non-zero 
values. T (N) and F (N) are continuously computed from the passive system, 
and TSET (N) and RATE (N) are controlling system characteristics supplied as 
initial constants. 
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Table 9 


Definition of Symbols Used in the Description of the Controlling System 


Symbol 

Vector 

Length 

Definition 

Dimension 

TSET (N) 

25 

"Set point" or reference point for receptors 
in compartment N 

°C 

ERROR (N) 

25 

Output from thermoreceptors in compartment N 

°C 

RATE (N) 

25 

Dynamic sensitivity of thermoreceptors in N 

h 

COLD (N) 

25 

Output from cold receptors in N 

°C 

WARM (N) 

25 

Output from warm receptors in N 

°c 

COLDS 


Integrated output from skin cold receptors 

°C 

WARMS 


Integrated output from, skin warm receptors 

°C 

SWEAT 


Total efferent sweat command 

Kcal . h" 1 

CHILL 


Total efferent shivering command 

Kcal . h" 1 

DILAT 


Total efferent skin vasodilation command 

1 . h _1 

STRIC 


Total efferent skin vasoconstriction command 

i 

N. D. 

SK3NR (I) 

6 

i 

Fraction of all skin receptors in Segment I 

N. D. 

SKINS (I) 

6 

i 

Fraction of sweating command applicable to 
skin of Segment I 

N. D. 

SK3NV (I) 

1 

6 

I 

Fraction of vasodilatation command applicable 
to skin of Segment I 

N. D. 

SKINC (I) 

i 

6 

Fraction of vasoconstriction command applicable 
to skin of Segment I 

N. D. 

MWORK (I) 

6 

Fraction of total work done by muscles in 
Segment I 

N. D. 

MCHIL (I) 

6 

Fraction of total shivering occurring in muscles 
of Segment I 

N. D. 
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Table 


1 0 


Definition of Symbols Used as Control Coefficients 


Symbol 

Definition 

Dimension 

csw 

Sweating from head core 

Kcal . h" 1 . °C” 1 

ssw 

! 

Sweating from skin 1 

Kcal . h" 1 . °C~ 1 

CDIL 

i 

Vasodilatation from head core 

1 . h -1 . “C" 1 

SDIL 

Vasodilatation from skin 

1 . h _1 . °C _1 

CCON 

Vasoconstriction from head core 

°c _1 

SCON 

Vasoconstriction from skin 

’c- 1 

CCHIL 

Shivering from head core 

Kcal . h" 1 . °C _1 

SCHIL 

Shivering from skin 

Kcal . h _1 . °C _1 

PSW 

Sweating from skin and head core 

Kcal . h" 1 . °C' 2 

PDIL 

Vasodilatation from skin and head core 

-1 -2 
1 1 . h . °C 

PC ON 

Vasoconstriction from skin and head core 

CM 

1 

o 

o 

PCHIL 

Shivering from skin and head core 

n -2 
1 °C 

BULL 

Factor determining temperature 

i 

°c -1 


sensitivity of sweat gland response 
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There can be conceptual controllers (1, 21) which are based on outputs which 
have a positive value on one side of the reference point, and have no output on 
the other side of the set point. This conceptual construction is supported by the 
obvious fact that there can be no negative firing rate, or negative sweating, 
negative blood flow or negative shivering. Non-linearity is thus introduced some- 
where in the chain from receptors to effectors, and it is convenient to have this 
non-linearity available at the receptor level. This is conveniently accomplished 
by testing ERROR (N) for its sign. If it is positive, we can assume warm receptors 
are active, and WARM (N) can be set to the value of ERROR (N) . If ERROR (N) is 
negative, we are dealing with a cold receptor output, and COLD (N) can be set to 
-ERROR (N). The resulting lists of WARM (N), COLD (N) and ERROR (N) then con- 
tain all thermal information which could be available to the thermoregulatory 
controller. There is very little information concerning the distribution of receptors 
in different compartments. This is of special importance in tissues which are wide- 
spread so that different compartments can be at widely different temperatures. The 
skin is a prime example. In the absence of any further information, the simplest 
approach would be the weight of the different compartments in different segments 
according to their mass. For the skin, Aschoff and Wever (15) have proposed rela- 
tive sensitivities which can be adopted although their factual basis is not clear. 
Bader and Macht (2 8) show a relatively higher sensitivity in the skin of head and 
chest than in the extremities, but there is reason to suspect that in their case the 
heat applied to the skin was sufficient to warm the central blood appreciably so 
that deeper or central thermoreceptors could have been involved in producing the 
results they observed. 
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There is little doubt that the trunk skin plays a more than proportionate role in 
skin thermoreception as evidenced by the difference in response to a uniform skin 
temperature and to an identical average skin temperature with divergent local tem- 
peratures. The tentative values for SKINR (I) shown in Table 11 have been obtained 
by multiplying the number of cold points for various segments given by Aschoff and 
Wever (15) with the corresponding surface areas and assigning to each area the 
appropriate fraction of the total sum of these products. 

Hertzman and Randall (18) surveyed the skin with respect to local blood flow. 
Although their measurements were photoelectric and consequently hard to quantitate, 
the relative proportions can probably be accepted. After rearrangement, the follow- 
ing fractions of the total blood flow are obtained for the basal and maximally dilated 
state respectively: head 0.288 and 0.132, trunk 0.265 and 0,322, arms 0.078 and 
0.095, hands 0. 10 and 0. 122, legs 0. 186 and 0.23, and feet 0.083 and 0. 10. 

The fractions for the maximally dilated state are taken as the distribution of 
vasodilatation over the skin, SKINV (I), in Table 11. Vasoconstriction is also not 
uniform over the skin area. In addition, the effects of countercurrent heat exchange 
will tend to exaggerate the effects of vasoconstriction in the distal extremities. 

On the other hand, countercurrent heat exchange is not so pronounced in the skin 
of the head and the trunk. Based on this, the assignments of SKINC (I) are made 
as shown in Table 11. 

The relative distribution of sweat secretion in the absence of information based 
on local secretion measurements is perhaps best approximated by the distribution 
of sweat glands. Randall (29) made such observations and the different skin areas 
were weighted as to the number of sweat glands, and the fractional distribution is 
given as SKINS (I) in Table 11. In addition, some weight was given to provisional 
results from thermograms presented by Stolwijk, et al. (30). 
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Total 1.8877 100.00 1.0000 1.000 1.000 

















The total warm receptor output from the skin WARMS can be obtained by 
* * 

summing SKINR (I) WARM (4 I) for the skin compartments of all six segments (I). 
The total cold receptor output from the skin can be integrated from a similar 
summation. 

It depends on the concept whether the signals from skin warm receptors are 
integrated without being offset by simultaneous cold receptor output from other 
skin areas. Available data is not adequate to clearly determine where integration 
of cold receptor signals and warm receptor signals takes place, so that the regu- 
lator in the controlling system must have either form available. Conceptual 
regulators which rely on linear addition of thermal signals from different sets of 
receptors (10, 20) use the ERROR (N) value in the controller equations, whereas 
conceptual models which use non-linear thermal signals either in a linear (10) 
or a non-linear expression ( 1 ) should use the WARM (N) or COLD (N) values 
in their controller equations. It will be clear that the weighted skin signal can 
be recovered in its linear form by subtraction of COLDS from WARMS. 

In the suggested controller equations which follow, we have assumed that 
the central receptors are located in the brain, so that WARM (1), COLD (l) or 
ERROR (1) constitutes the thermal signal from the interior. There are reasons to 
believe that central receptors occur elsewhere (25, 28) but not enough is known 
to take them separately into account, and central temperatures will not normally 
be found to move far apart, so that lumping central receptors into the brain com- 
partment will only result in a second order error for most circumstances. 

In general, the controller equations produce signals which drive the effector 
part of the regulator. These drives are not necessarily identical to effector action 
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since they may be modified at the effector site by mechanisms responding to local 
conditions. The controller equations all have a term consisting of the product of 
a control coefficient and a central temperature signal, a term consisting of the 
product of a control coefficient and an integrated skin temperature signal and a 
third term consisting of the product of a control coefficient, a central temperature 
signal, and a skin temperature signal. Controller equations are then as follows: 

SWEAT = CSW* ERROR (1) + SSW* (WARMS -COLDS) + PSW* WARM (1)* WARMS 

DILAT = CDIL* ERROR (1) + SDIL* (WARMS - COLDS) + PDIL* WARM (1)* WARMS 

CHILL = -CCHIL* ERROR (l) + SCHIL* (COLDS-WARMS) + PCHIL* COLD (1)* COLDS 

STRIC = -CCON* ERROR (1) + SCON (COLDS - WARMS)+ PCON* COLD (1)* COLDS 

Since these expressions can become negative under certain circumstances, it 
is necessary to make sure that any negative values of SWEAT, DILAT, CHILL or STRIC 
are set to zero. 

The third part of the controller distributes the effector commands to the appropriate 
compartment and, if required and appropriate, modifies them according to local 
conditions. 

In order to facilitate the description of this distribution, the following con- 
vention will be used: for segment I, the core layer will be indexed as N (equal 
to 4 1-3), which makes the concentric muscle layer be indexed as (N + 1); the 
subcutaneous fat layer and the skin layer for segment I then are indicated by the 
subscripts (N + 2) and (N + 3) respectively. The metabolic rates for all compartments 
then become as follows: 

Q (N) = QB (N) 
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This expression states that the basal metabolic rate assigned to the cores of 
the various segments is not thought to vary under the relatively short term con- 
ditions for which the model will be used. 

Q (N + 1) = QB (N + 1) + WORKM (I)* (WORK!) + CHILM (I)* CHILL 
The metabolic heat production in the muscle compartment is the sum of the 
basal rate of heat production QB (N + 1) , and the sum of the heat production 
rates assigned due to muscular work done and any shivering thermogenesis. 
WORKI is the total heat production as one would determine from the oxygen up- 
take minus the basal metabolic rate and the caloric equivalent of the external 
work performed. WORKM (I) and CHILM (I) are estimates of the distribution of 
extra heat production over the various muscle compartments. Th estimates used 
are given in Table 12. 

Table 12 


Estimates of Distribution of Heat Production in Muscle Compartments 



McOTa.. 

WORKM (I) 

CHILM (I) 

Head 

2.323 

— 

0.02 

Trunk 

54.790 

0.30 

0.85 

Arms 

1 0.525 

0.08 

0.05 

Hands 

0.233 

0.01 

0.00 

Legs 

31.897 

0.60 

0.07 

Feet 

0.233 

0.01 

0.00 
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In the absence of data, the hec '- from shivering wac distributed according to the 
mass in the various muscle compartments, whereas the values for WORKM (I) are 
based on estimates for bicycle exercise. Values for WORKM (I) should be reviewed 
for other types of exercise. 

For the subcutaneous fat and the skin compartments, it is assumed that the 
basal heat production does not change in the conditions to be evaluated by the 
model, so that: 

Q (N + 2) = QB (N + 2) 

Q (N + 3) = QB (N + 3) 

The convective heat transfer by blood flow plays a very important role in the 
thermal responses to internal and external stresses. Certain simplifications are 
implicit in the expressions which follow. Some of these simplifications are due 
to a scarcity of applicable data, some of them are justified by the fact that although 
inaccuracies in mass flow are introduced by the simplifications, these inaccuracies 
have only a very small effect on heat transfer. For example, it is assumed that 
blood flow to the core compartment remains at the basal value: 

BF (N) = BFB (N) 

This ignores the fact that blood flow to the trunk core in fact can be reduced 
during exercise stress (16, 31, 32), but, in fact, the blood flow to these compart- 
ments is considerably in excess of that required to supply oxygen to these compart- 
ments, so that only small thermal gradients can ever exist between the head core, 
trunk core and the central blood. Small reductions in core blood flow will then re- 
sult in slight increases in these small gradients, and convective heat transfer is not 
materially affected. The muscle compartments, on the other hand, can have widely 
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varying metabolic rates, and substantial variations in blood flow. Regional blood 
flow measurements overworking quadriceps muscle (34) and measurements of ar- 
teriovenous oxygen differences across the working quadriceps muscle (35) can be 
added to measurements of local muscle temperature (36) during exercise to arrive 
at an insight of muscle blood flow during exercise. Under conditions where 
heat loss in exercise is mostly from evaporation of sweat, the temperature of the 
working muscle is about 1.0 °C above that of the arterial blood supplying oxygen 
(35). If venous blood leaves at muscle temperature, then every liter of blood 
carries off 0.95 Kcal. Each liter of blood contains about 200 ml 0 2 which, if 
completely scavenged, could produce about 1.0 Kcal of heat. It is thus reason- 
able to assume that at least 1 1 of blood is required to produce 1 Kcal of heat. 

Asa result, the following equation describes the blood flow to the muscle compart- 
ments . 

BF (N + 1) = BFB (N + 1) + Q (N + 1) - QB (N + 1) 

Blood flow to the subcutaneous fat is not very high in basal value, and it is 
assumed that it is not effectively changed as a result of thermoregulatory adjustments: 

BF (N + 2) = BFB (N + 2) 

Skin blood flow is highly dependent on the thermoregulatory controller. The 
basal blood flow at thermal neutrality can be reduced to very low values through 
vasoconstriction, and increased substantially through vasodilatation; in addition, 
there is an effect of the local temperature which has a substantial modifying influ- 
ence on resistance of cutaneous veins,' as demonstrated by Webb-Peploe and 
Shepherd (3 6). These investigators showed that the local temperature of the 
saphenous vein had a very strong effect on the response of the vein to a given level 
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of constricting stimulation. The peripheral resistance caused by a given level of 
sympathetic tone would be reduced by high local temperatures and increased by low 
local temperatures. In many cases, the resistance would about double for a 10°C 
drop in local temperature, and for a given central signal. The tentative expression 
for skin blood flow then becomes: 

ic "k * * * 

BF (N) = ( (BFB (N) + SKINV (I) DILAT) / (1. + SKINC (I) STRIC) ) 2. (ERROR (N)/10.) 

In this expression, the weighted DILAT drive is added to the basal blood flow. 

The weighted constrictive tone operates through a resistance and thus is entered as 
a divisor. The local skin temperature effect then acts on the total flow drive multi- 
plying it by unity in neutral conditions, by less than unity at skin temperatures below 
normal and by more than one at skin temperatures above the neutral value. 

It could be that the local effect on cutaneous veins described by Webb-Peploe 
and Shepherd does not affect total skin blood flow very much, but reduces heat loss 
by shifting venous return to deeper veins. This shift would activate an effective 
counter-current heat exchange between arteries and veins, which could not be dis- 
tinguished from a reduction in skin blood flow as far as convective heat exchange 
is concerned. 

Evaporative heat loss occurs from the trunk through respiratory water loss, and 
through evaporation from the skin. Respiratory water loss is a function of the water 
vapor pressure in the inspired air (PAIR) and the ventilation volume. The ventilation 
volume is closely related to the metabolic rate, and ventilatory evaporative heat loss 
can be approximated over a wide range of meta-bolic rates by: 

E (5) = (74.4 + WORKI)* 0.023* (44.0 - PAIR) 
in which 44.0 is the vapor pressure in expired air. 
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For all other core compartments 

E (N) = EB (N) 

and since muscle and subcutaneous fat compartments have no evaporative heat 
loss: 

E (N + 1) = EB (N + 1) 

E (N + 2) = EB (N + 2) 

The different skin compartments receive a sweating drive from the central 
controller. The local response in a given compartment depends on the skin 
area and on the number of sweat glands present, expressed by SKINS (I). In 
addition, the local skin temperature modifies the response to a given sweat 
command as described by Bullard and associates (24). 

On repeating these observations in our laboratory, we confirmed the local 
temperature effect, but found it to have a somewhat lower as shown in 
Figure 2. As a result, the of the local temperature effect which Bullard, et al. 

report as about 8 is used here at a value of 2. The evaporative heat loss of the 
skin of each segment then becomes: 

E (N + 3) = (EB (N'+ 3) + SKINS <I>* SWEAT)* 2.** ( ,T (N * 3) ‘ TSET (N + 3) )/BULU 
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SWEATING RATE mg/cm 2 min 



Figure 2. Dependence of local sweat rate on the 
skin of the thigh, on average skin tem- 
perature, and local skin temperature. 
Internal temperatures were constant during 
these measurements. 
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GOO 


Annotated Listing of FORTRAN Statements Implementing the Model 

A flow sheet of the program for simulation of body temperature regulation in man 
is given in Figure 3, and the listing will follow by sections as identified in the flow 
diagram. 

DIMENSION T(25), TSET (25), RATE (25), C(2S), QB(24), EB(24), BFB (24) 
DIMENSION TC(24), S(6), SHNR(6), SKINS (6), SKINV(6), SKINC (6), WORKM (6) 
DIMENSION CHILM (6), HR (6), HC(6), P (10), F (25), H(6), WARM (25), COLD(2S) 
DIMENSION HF (25), PSKIN (6) 

DIMENSION ERROR (25), Q (24). E (24), BF (24), EMAX(6), BC (24), TD(24) 

READ CONSTANTS FOR THE CONTROLLED SYSTEM 

100 FORMAT (14F5. 2) 

101 CONTINUE 


READ (2, 100) 

C 

READ (2, 100) 

QB 

READ (2, 100) 

EB 

READ (2, 100) 

BFB 

READ (2, 100) 

TC 

READ (2, 100) 

S 

READ (2, 100) 

HR 

READ (2, 100) 

HC 

READ (2, 100) 

P 

SA = 0.0 
DO 110 K = 

1, 6 


110 SA = SA + S (K) 
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READ CONSTANTS 


READ 


READ 


READ 

CONTROLLED 

> 

CONTROLLER 

— > 

INITIAL 

B5 

EXPERIMENTAL 

SYSTEM 


CONSTANTS 


CONDITIONS 

i 

CONDITIONS 
1 


PREPARE 

OUTPUT 


t 

OUT 

PUT 


YES 



CALCULATE 

NEW 

TEMPERATURES 


ESTABLISH 

THERMORECEPTOR 

OUTPUT 


INTEGRATE 

PERIPHERAL 

AFFERENTS 


DETERMINE 

EFFERENT 

OUTFLOW 


DETERMINE 


CALCULATE 


ASSIGN 

INTEGRATION 

4— 

HEAT FLOWS 

< 

EFFECTOR 

STEP 




OUTPUT 


Fig. 3. Flow diagram for simulation of thermoregulation 

in man 
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The first section reads in the constants which define the controlled system as 

defined in the tables before. The values for P are taken from the steam table and 

consist of the saturated water vapor pressures at 5°C intervals from 5°C to 50°C. 

2 

Statement 110 calculates the total surface area SA in m for use in the output 
preparation phase. 

The second phase reads in the controller constants beginning with the set point 
temperatures. As set point temperatures, the values of steady state equilibrium 
temperatures are used which are reached under the following conditions: air tem- 
perature 29.4S°C, air velocity 0.1 m/sec, relative humidity 0.3, no work, and all 
controller coefficients set to 0.0, with the exception of BULL which is set to 10.0. 
C 

C READ CONSTANTS FOR THE CONTROLLER 
C 


READ (2, 100) 

TSET 

READ (2, 100) 

RATE 

READ (2, 100) 

SKINR 

READ (2, 100) 

SKINS 

READ (2, 100) 

SKENV 

READ (2, 100) 

SKINC 

READ (2, 100) 

WORKM 

READ (2, 100) 

CHILM 

READ (2, 100)CSW, SSW, PSW, CDIL, SDIL, PDIL, CCON, SCON, PCON, 

CCH1L, SCH1L 

XPCHIL, BULL 
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The values for SKINR, SKINS, SKINV, SKINC, WORKM and CH1LM have been de- 
fined before in preceding Tables 11 and 12. 

C 

C READ INITIAL CONDITIONS 
C 

READ (2, 100) T 
TIME = 0.0 
ITIME = 0 
JTIME = 0 
DO 102 N = 1, 25 
F (N) =0.0 
102 CONTINUE 

The initial conditions for the next phase consist of the initial values for the 
temperatures in all compartments. These values are the same as those for TSET 
and were obtained as described there. The values are given in Table 13. 

C 

C READ EXPERIMENTAL CONDITIONS 
C READ (2,299) TAIR, V, RH, WORK, INT 

299 FORMAT (4F5.0, 12) 

IF (WORK - 74.4) 104, 104, 105 

104 WORKI = 0.0 
GO TO 106 

105 WORKI = (WORK - 74.4) * 0.78 

106 CONTINUE 

DO 202 1 = 1,6 

H(I) = (HR(I) + HC (I) * (V/O.l) **0.5) * S (I) 

202 CONTINUE 
I = TAIR/5 

PAIR = RH*(P(I) + (P (I + 1) - P (I) )*(TAIR - 5 * I)/5.0) 
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Table 


1 3 

Set Point Values and Initial Condition Temperatures 


Segment 

Compartment 

N 

Temperature °C 

Head 

Core 

1 

36.96 


Muscle 

2 

35.07 


Fat 

3 

34.81 


Skin 

4 

34.58 

Trunk 

Core 

5 

36.89 


Muscle 

6 

36.28 


Fat 

7 

34.53 


Skin 

8 

33.62 

Arms 

Core 

9 

35.53 


Muscle 

10 

34.12 


Fat 

11 

33.59 


Skin 

12 

33.25 

Hands 

Core 

13 

35.41 


Muscle 

14 

35.38 


Fat 

15 

35.30 


Skin 

16 

35.22 

Legs 

Core 

17 

35.81 


Muscle 

18 

35.30 


Fat 

19 

35.31 


Skin 

20 

34 .10 

Feet 

Core 

21 

35.14 


Muscle 

22 ! 

35.03 


Fat 

23 

35.11 


Skin 

24 1 

1 

35.04 

Central Blood 


25 

36.71 

























The next phase reads in the experimental conditions, air temperature, air velocity 
in m. sec *, the relative humidity as a fraction rather than a percentage, the work 
rate in total heat production in Kcal.h * as normally computed from oxygen uptake, 
and the interval in minutes between desired outputs. Next, WORKI is computed by 
subtracting the basal metabolic rate, and by subtracting the external work, so that 
WORKI represents the total extra heat production in the working muscle. Subse- 
quently, H(I) is calculated in Kcal.h * . °C * for each segment. The last two 
statements look up the partial pressure of water vapor in the environment by inter- 
polation from the steam table. 

C 

C ESTABLISH THERMORECEPTOR OUTPUT 
C 

301 CONTINUE 

DO 302 N = 1, 25 
WARM (N) = 0.0 

COLD (N) = 0.0 

IF (F(N) ) 310, 311, 311 
311 F (N) = 0.0 

310 CONTINUE 

ERROR (N) = T (N) - TSET (N) + RATE (N) * F (N) 

if (Error (n) ) 303 , 302 , 304 

303 COLD (N) = ERROR (N) 

GO TO 302 

304 WARM (N) = ERROR (N) 

3 02 CONTINUE 
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In the next phase, the actual instantaneous temperatures in all compartments are 
compared with the "set point" temperatures, taking into account the rate of change 
in those cases where this is appropriate. The direction of rate of change is first 
determined, and positive rates of change are set to zero. This was done because 
we have not been able to demonstrate specific effects of a positive rate of change, 
although a direct effect of negative rate of change of skin temperature on sweating 
has been demonstrated. Such effects had been implied earlier (3 7) and were studied 
quantitatively in our laboratory by Nadel, Bullard and Stolwijk (3 8). In these exper- 
iments it was shown that local sweating rate as measured by a sweat capsule and 
resistance hygrometry was linearly dependent on average skin temperature, as long 
as local skin temperature and internal temperatures were held constant. Radiant 
energy was used to bring about rapid changes in average skin temperature. Figure 
4 shows such linear dependence on average skin temperature at various local skin 
temperature, in a steady state condition. When rapid changes were brought about, 
it was found that during skin warming the instantaneous mean skin temperature was 
a good predictor of instantaneous local sweat rate, but during cooling the sweat rate 
at any instant was lower than that predicted by instantaneous mean skin temperature. 
It was found that this deviation was proportional to the negative rate of change in 
average skin temperature as shown in Figure 5. 
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Figure 4. 


Figure 5. 



Linearity of local sweating response with average skin 
temperature, at constant internal temperature. 


DEVIATION FROM 
LINEAR SWEATING RESPONSE 

mg / cm* min 

♦ 0.20 -j~ 


dT, 
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C per min 




*0.5 * *IO 


Deviation from sweating response linear with average skin 
temperature, as a function of rate of fall of mean skin temperature. 
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This phase results in lists of temperature deviations from normal, in linear 
terms (ERROR (N) ) and in non-linear terms WARM (N) and COLD (N) which are 
zero for values below and above the set point respectively, 

C 

C INTEGRATE PERIPHERAL AFFERENTS 
C 

WARMS = 0.0 
COLDS = 0.0 
DO 305 I = 1, 6 
K = 4*1 

WARMS = WARMS + WARM (K) * SKINR (I) 

COLDS = COLDS + COLD (K) * SKINR (I) 

305 CONTINUE 

Since only the skin afferents are used in the present model, these signals are 
the only ones to be weighted for integration. Weighting is done according to the 
receptor density estimates SKINR (I). The use of WARM (K) and COLD (K) does not 
necessarily cause the loss of linear information. A linear form of the weighted 
value of the skin error signal can be recovered by the use of WARMS - COLDS. 

C 

C DETERMINE EFFERENT OUTFLOW 
C 

SWEAT = CSW*ERROR(l) +SSW* (WARMS -C OLDS) +PSW* WARM (1) * WARMS 
DILAT = CDIL* ERROR (1) + SDIL* (WARMS - COLDS) + PDIL* WARM (1) * WARMS 
STRIC = -CCON* ERROR (1) -SCON* (WARMS - COLDS) + PCON * COLD (1)* COLDS 
CHILL = -CCHIL* ERROR (1) -SCHIL* (WARMS-COLDS) + PCHIL * COLD (1) * COLDS 
IF (SWEAT) 309, 312, 312 
309 SWEAT = 0.0 
312 IF (DILAT) 313, 314, 314 
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313 DILAT 


0.0 


314 IF (STRIC) 315, 316, 316 

315 STRIC = 0.0 

316 IF (CHILL) 317, 318, 318 

317 CHILL =0.0 

318 CONTINUE 

The next section calculates the control commands going to the periphery. 

Since the ERROR (N) terms can become negative as well as positive, each of the 
commands is protected against becoming negative. 

ASSIGN EFFECTOR OUTPUT TO THE COMPARTMENTS 

400 CONTINUE 

DO 401 I = 1, 6 
N = 4*1-3 
Q (N) = QB (N) 

BF (N) = BFB (N) 

E (N) = EB (N) 

Q (N + 1) = QB (N + 1) + WORKM (I) * WORK3 + CHILM (I) * CHILL 
E (N + 1) = 0.0 

BF <N + 1) = BFB (N + 1) + Q (N + 1) - QB (N + 1) 

Q (N + 2) = QB (N + 2) 

E (N + 2) = 0.0 
BF (N + 2) = BFB (N + 2) 

Q (N + 3) = QB (N + 3) 

E(N + 3) = (EB(N+3) + SKINS (I)* SWEAT) *2.0** (ERROR (N + 3) /BULL) 

BF (N + 3) = ((BFB(N + 3)+SKINV(I) * DILAT)/(1 . 0 + SKINC (I)* STRIC) )* 2 . 0** 
X(ERROR (N + 3)/ 10.0) 
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K = T (N + 3)/5 

PS KIN (I) = P (K) + (P (K + 1) -P(K) ) * (T (N + 3) -5*K) / 5 . 0 
EMAX (I) = (PS KEN (I) -PAIR) *2 . 2 * (H (I) -HR (I) * S (I) ) 

IF (E (N + 3) -EMAX (I)) 403, 403, 402 

402 E (N + 3) = EMAX (I) 

403 CONTINUE 

401 CONTINUE 

E (5) = (74.4 + WORM) * 0.0023 * (44.0 -PAIR) 

This phase distributes the control commands to the various segments and within 
the segments to compartments in the form of the calculated effect of heat flow in 
Kcal . h 1 for Q(N) and E(N), and of conductance in Kcal . h 1 . °C 1 for BF (N) . 

The core compartments (N) and the fat compartments (N + 3) remain at basal values. 

The muscle compartments receive their assigned portion of the total internal heat 
production due to work, and of that due to shivering. For each Kcal of heat produced, 
the muscle requires at least one liter of blood, as seen in BF (N + 1). 

Skin metabolic rate remains at basal values, but skin blood flow and evaporative 
heat loss are determined by command signals, as modified by the local temperature. 

The water vapor pressure at the skin surface temperature is read and interpolated 
from the steam table and the maximum evaporative heat loss from a completely wet 
skin EMAX (I) is computed, based on existing conditions . If the computed evaporation 
from the skin E (N + 3) exceeds the maximum possible EMAX (I), then the former is 
set equal to EMAX (I). The last statement calculates the evaporative heat loss from 
the respiratory tract with the aid of total metabolic rate and the water vapor pressure 
in the environment. 
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c 

C CALCULATE HEAT FLOWS 
C 

DO 500 K = 1, 24 

BC (K) = BF (K) * (T (K) -T (25) ) 

TD (K) = TC (K) * (T (K) -T (K + 1) ) 

500 CONTINUE 

DO 501 I = 1, 6 
K = 4*1-3 

HF (K) = Q (K) -E (K) -BC (K) -TD (K) 

HF (K + 1) = Q (K + 1) -BC (K + 1) + TD (K) -TD (K + 1) 

HF (K + 2) = Q (K + 2) -BC (K + 2) + TD (K + 1) -TD (K + 2) 

HF (K + 3) = Q (K + 3) -BC (K + 3) -E (K + 3) + TD (K + 2) -H(I) * (T(K+3)-TAIR) 

501 CONTINUE 
HF (25) = 0. 

DO 502 K = 1, 24 
HF (25) = HF (25) + BC (K) 

502 CONTINUE 

The next section computes the net heat flow rates into or out of each of the com- 
partments. Initially BC (K), the convective heat flow rate, and TD (K), the conductive 
heat flow rate for each compartment to the next compartment are calculated. After 
this, all the heat flow rate components are added up for each compartment. The only 
exception is the central blood compartment which receives the sum of all convective 
heat flows. 
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c 

C DETERMINE OPTIMUM INTEGRATION STEP 
C 

DT = 0.016666666667 
DO 600 K = 1, 25 
F (K) = HF (K)/C 00 
U = ABS (F (K) ) 

IF (U * DT -0. 1) 600, 600, 601 

601 DT = 0.1/U 

600 CONTINUE 

The optimum integration time step is determined in the next phase. Initially 
the time increment for numerical integration is set to 0.01667 hours or 1 minute 
this being considered the shortest plausible print-out interval. Based on this 
increment, the temperature steps in each compartment are calculated and if any 
exceed 0.1 °C then the time increment DT is reduced so that the maximum tem- 
perature change in any compartment in each iteration is kept to 0. 1°C or less. 

C 

C CALCULATE NEW TEMPERATURES 
C 

DO 700 K = 1, 25 
T (K) = T (K) + F (K) * DT 

700 CONTINUE 

TIME = TIME + DT 
LTIME = 60. * TIME 

IF (LTIME -INT -ITIME) 301, 701, 701 

701 CONTINUE 
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Based on the time increment thus arrived at, the next section calculates 
the new temperatures and increments the clock with the value of the integration 
increment DT. This followed by a test to see whether output is desired. If such 
output is not desired, the program returns to the phase where thermoreceptor output 
is re-computed. If output is desired, the program proceeds to the next phase. 

PREPARE FOR OUTPUT 
EWET = 0.0 
DO 450 1 = 1,6 

EWET = EWET + (E (4 * I)/EMAX (I) ) * (S (I)/SA) 

450 CONTINUE 

PWET = 100. * EWET 

I TIME = I TIME + INT 
CO = 0.0 
HP = 0.0 
EV = 0.0 
TS = 0.0 
TB = 0.0 
HFLOW = 0.0 
SBF =0.0 
DO 800 N = 1, 24 
CO = CO + BF (N)/ 60.0 
HP = HP + Q (N) 

EV = EV + E (N) 

800 CONTINUE 

DO 802 I = 1, 6 
SBF = SBF + BF (4 * I)/ 60.0 
TS = TS + T (4 * I) * C (4 * D/3.35 


46 



802 


CONTINUE 


DO 801 N = 1, 25 

TB = TB + T (N) * C (N)/59.17 

HFLOW = HFLOW + HF (N) 

801 CONTINUE 

EV = EV/SA 
HP - HP/SA 
HFLOW = HFLOW/ SA 

COND = (HP -E(5)/SA -HFLOW) / (T (25) -TS) 

In this phase, the value of variables which refer to compartments or segments 

are combined to yield values in the dimensions usually used for reporting thermo- 

physiological experiments. The first few statements calculate the fractional wetted 

area of the skin according to Gagge (13) from the area-weighted ratios of the actual 

evaporation rate and the maximum possible evaporation rate. The output time clock 

is incremented with one print-out interval time. The cardiac output, heat production 

and evaporative heat loss are obtained by summing of blood flow, metabolic heat 

production rate and evaporative heat loss over all compartments. Skin blood flow, SBF, 

and mean skin temperature, TS, are calculated by summing of segmental skin blood flows 

and skin temperatures. Similarly, mean body temperature is obtained by averaging 

of all compartmental temperatures weighted with their thermal capacitance. Net 

rate of heat storage for the whole body is obtained by summing all net heat flow over 

-2 -1 

all compartments. EV, HP and heat flow are reduced to Kcal . m . h by dividing 
them by the total body surface area SA. The core-to-skin thermal conductance, COND, 
which is often derived in physiological experiments is derived in the conventional 
manner. 
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The next phase contains the output statements which are not a material part of 
the program and which would normally depend on the purpose for which the simula- 
tion is intended. 

BEGINNING OF OUTPUT SECTION 
CALL DATSW (0, K) 

GO TO (951, 950), K 
951 CONTINUE 

IF (ITIME -INT) 909, 909, 911 

909 PAUSE 

910 WRITE (1, 912) 

912 FORMAT ('TIME S M EV TB TS TH TO TR TM 1 SBF CO COND PWET') 
NN = 0 

911 IF (NN-42) 913, 913, 914 

913 WRITE 915) ITIME, HFLOW, HP, EV, TB, TS, T(l), T(25), T(5), T(18), SBF, 
CO, 1 COND, PWET 

915 FORMAT (13, 3F7.1, 8F6.2, 2F6.1) 

NN = NN + 1 
GO TO 1100 

914 WRITE (1, 916) 

916 FORMAT (22(/) ) 

GO TO 910 

950 CONTINUE 
1100 JTIME = JTIME+INT 
CALL DATSW (1, K) 

GO TO (917, 1102), K 
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917 CONTINUE 


ITA = TAIR * 10. 

ISKSW = 1.16 * (EV-9.5-0.08 * WORK/SA) 

IHP = HP 
ITS = TS * 100. 

ITR = T (5) * 100. 

ITO = T (25) * 100. 

ITM = T (18) * 100. 

IEV = EV 

ICOND = COND *1.16 
IRM = HP * 1. 16 
ITB = TB * 100. 

WRITE (2, 9901) ITIME, ITA, IHP, ITS, ITR, ITO, ITM, IEV, ISKSW, ICOND, 
IRM, XITB, SBF, CO 

9901 FORMAT (313, 414, 3X, 13, 3X, 313, 4X, 14, F4.2, F5.2) 

1102 IF (JTIME-30) 301, 1101, 1101 
1101 JTIME = 0 
PAUSE 

IF (ITIME- 2 10) 102, 101, 101 
9990 CALL LEAVE 
901 CALL EXIT 
END 
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Control Coefficients 


As will be evident from simulations shown later, it is not yet possible 
to account for all experimental results with a single set of control coefficients. 

This section is concerned with the description of current controller concepts and 
values of control coefficients and the experimental justification and validation. 

It will be seen that there is a tendency to obtain values from a variety of steady 
state responses with backgrounds as diverse as possible, and to attempt to 
validate them under dynamic conditions. The dynamic conditions do not only test 
the validity of concept and quantitative character of the controller, but they also 
challenge the validity of the controlled system as described here. The relatively 
small number of nodes in the finite difference method used here does impose a 
limit on the time resolution which can usefully be obtained, especially in con- 
ditions where large gradients exist, or where vasoconstriction reduces convective 
heat transfer. 

Control of Sweating 

We have available in the results obtained in this laboratory in resting ex- 
periments and exercise experiments the basic information to describe the control 
of sweating. Figure 6 shows a composite plot of more or less steady state rates 
of sweating as dependent on internal temperature and skin temperature, taken 
from (10, 11, 30). The broken line demarcates the area outside which, for physical 
or physiological reasons, it would be hard or impossible to produce a steady state. 

It is clear that the different lines of equal mean skin temperatures have 
different slopes as well as intercepts. From Figure 6 and Figure 4, we can then 
draw the following conclusions: 
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1. At constant skin temperature, sweating is proportional 
to internal temperature (Figure 6). 

2. At constant internal temperature, sweating is proportional to 
mean skin temperature (Figure 4) . 

3. At a given combination of internal temperature and mean skin 
temperature, the local sweating is dependent on local skin 
temperature with a of about 2 (Figure 4). 

It is clear then that Figure 6 has included in it this local effect. By making 
the arbitrary assumption that the mean skin temperature at thermal neutrality 
is the "normal" one for referencing the peripheral Q effect, we can eliminate 
the local effect from Figure 6 and this produces Figure 7 which then represents 
the output from the central controller. The expression of the controller output 
in Figure 7 then becomes: 

SWEAT = 320. * ERROR (I) + 29. * (WARMS - COLDS) so that 

CSW - 320.0 SSW = 29.0 and PSW =0.0 become the controller coefficients 

for the central controllers. From Figure 6, it can be derived that the local skin 

temperature effect has a Q 10 of 2 so that the local coefficient BULL = 10.0, and 

** 

the local multiplier becomes: 2 . 0 (ERROR (N + 3) / 10. 0) . 

This concept of the control of sweating thus contains the additive as well 
as multiplicative aspects of integration of central and peripheral signals, and, 
moreover, does not require for steady state conditions a non-thermal factor 
which has been postulated. 
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Figure 6. Composite plot of more or less steady state rates of sweating 
as dependent on internal temperatures . 



Control of Vasodilatation 


Measurements of peripheral vasodilatation is not directly possible as is 

the case with sweating. The most widely used measurement is a derived one 

usually indicated by the symbol K which is the effective thermal conductance 

-2 -1 -1 

from core to skin surface, expressed in Kcal . m . h . °C ; K is cal- 
culated from: 

COND = (HP - E (5) - HFLOW) / ( (T (5) - TS)/SA) 

It is implicit in the definition that the thermal conductance can only be calculated 
during conditions where no gradients are changing; i.e., in quasi-steady states. 
Such values have more inherent fluctuation than direct measurements so that the 
plot assembling such values from a number of experimental conditions from (10, 

11, 30) is not as orderly as that for sweating. The plot is given in Figure 8. 

The similarity of this plot to that of Figure 6 is probably not altogether coinci- 
dental, in view of the similarity between the Bullard, et al. (24) findings for the 
action of peripheral temperature on sweat gland response, and the findings of 
Webb - Peploe and Shepherd (36) on the action of peripheral temperature directly 
on the responsiveness of cutaneous veins to sympathetic tone. 

If we may make the not yet completely proven assumption that the similarity 
is more than suggestive, it is useful to derive control equations for the central 
vasodilatory drive in a manner similar to that in the case of sweating. Figure 10 
then shows the value of K for different values of central temperature and mean 
skin temperature, with the local peripheral effect removed, again normalizing at 

the skin temperature at thermal neutrality. It is then found that the slope of the 

-2 -1 -1 -2 -1 
lines is 62 1 . m . h . °C if we equate an increase of 1 Kcal . m . h . °C 

- 1 — 9 

in the conductance with at least 1 1 . h . m increase in skin blood flow. It 
then also follows that 1 °C change in mean skin temperature at constant internal 
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temperature results in about 4 1 , m . h change in skin blood flow. The 
controller command then becomes: 

DILAT = 117, ERROR (I) + 7.5 * (WARMS - COLDS) 
so that CDIL = 117.0, SDIL = 7.5 and PDIL = 0.0 

It also becomes clear that the local effect causes a doubling or halving of 
local blood flow for a 6°C change in local temperature, so that the local mul- 
tiplier becomes: 

2.0 (ERROR (N + 3) / 6. 

The value for the central control coefficient CDIL of 117.0 is very close to 
one proposed in earlier reports and also compares with a value of about 90 
proposed by Benzinger, et al. (39). 

Control of Shivering 

Shivering, although easily measured, is among the elusive physiological 
responses since steady states are difficult to obtain. However, a large number of 
measurements all point to the fact that both skin temperature and internal tempera- 
ture must be below a certain threshold. Such measurements were made by Benzinger, 
et al. (39) and by Nadel (40). Probably as a result of differences in methods, 
these authors find slightly different values for the threshold for the central and 
mean skin temperatures. These authors, as do many others, report responses 
which can consistently be predicted by a multiplicative expression. Measurements 
in this laboratory have failed to give responses which were consistent enough to use 
for quantitative extrapolation. One very consistent difference was found in that 
for a given cold exposure women shivered earlier and to a much higher degree 
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Figure 8. Plot of a number of experimental conditions. 
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than men. Women also showed a more labile vasomotor response in heat and 
cold and had a consistently lower CSW value than men. The latter difference 
resulted in women reaching higher internal body temperatures during resting 
heat exposure. A detailed report on these experiments is in preparation. 

Benzinger, et al. (39 ) report results which would suggest that in the expression 
for shivering control the following value should be used: 

CHILL = 21.0 * COLD (I) * COLDS 
which would result in the following control coefficients: 

CCHIL =0.0, SCHIL =0.0 and PCHIL = 2 1 . 0 . 

Nadel’s results would indicate a higher value for PCHIL in the range of 40.0, but 
his measurements were short lasting and should, therefore, be approached with 
more caution. 

Control of Vasoconstriction 

The type of vasoconstriction implied here is largely limited to the special 
peripheral vascular beds in the hands and feet which normally have a high non- 
nutritive blood flow which is highly sensitive to both internal temperature and 
mean skin temperature. In addition, there is a high sensitivity to psychogenic dis- 
turbances. As the mean skin temperature decreases below normal, or as the central 
temperature decreases below normal, constrictor tone in hands and feet increase. 
Since this represents an increase in resistance, the calculated local blood flow is 
reduced by dividing it by a value STRIC. In the absence of clear data, the following 
expression is used to define STRIC: 

STRIC = -5.0* (ERROR (I) + WARMS - COLDS). 
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Performance of Model of Thermoregulation in Man at Rest 

It is easily recognized that the model will reproduce the steady state thermo- 
regulatory responses with acceptable accuracy since the control coefficients are 
essentially based on steady state relationships. The best method of validation 
then is to reproduce a sudden change in conditions with the resultant dynamic 
thermoregulatory response. In order to do this, two specific experiments were 
simulated which were also evaluated with an earlier model (1). Young male sub- 
jects were quickly transferred between two ambient temperatures; the air and 
wall temperatures were equal, and the air velocities were minimal (i.e., 0.1 m/sec 
or less). Average skin temperature, rectal temperature, tympanic membrane tem- 
perature, rate of evaporative weight loss and metabolic rate were determined con- 
tinuously. The experiments have been described in more detail before (10, 11). 

Figure 10 shows an experiment in which the subject spent 30 minutes at a thermally 
neutral temperature, followed by 120 minutes at an ambient temperature of 48°C, 
and an additional 60 minutes of recovery in a 30°C environment. In each case, 
the broken line represents the response predicted by the model, and the solid line 
shows the experimental result, in measurements of rectal temperature, mean skin 

temperature and sweat rate. The metabolic rate, during the whole experiment, stayed 

-2 -1 

in the range of 38 - 42 Kcal . m . h The results predicted by the model are given in 
more detail in Table 14. The headings in Table 14 and other tables following after Table 

14 are abbreviations for the following variables: TIME, in minutes; S: rate of heat 

-2 -1 -2 -1 
storage in Kcal . m . h ; M: metabolic heat production in Kcal . m . h ; 

-2 -1 

EV: total evaporative heat loss in Kcal . m . h ; TB: mean weighted body 
temperature in °C; TS: mean weighted skin temperature in °C; TH: temperature of 
the head core representing hypothalamic temperature in °C; TO: temperature of 
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the central blood compartment, representing esophageal temperature in °C; TR: 

trunk core temperature representing rectal temperature in °C; TM: temperature of 

the muscle compartment in the leg in °C; SBF: total skin blood flow in 1 . min 1 ; 

CO: estimated cardiac output in 1 . min COND: equivalent thermal conductance 

-2 -1 -1 

from core to skin, Kcal . m , h . °C ; PWET: percentage of the skin surface 
as wetted area (13, 14) for the sweat rate and environmental conditions given. 

In a different type of experiment, the subjects spent 60 minutes at an en- 
vironmental temperature of 43°C followed by 120 minutes at 17°C and a final 30 
minutes at 43°C. Data and simulation of this experiment are given in Figure 11 

and Table 15. In the experiment, the metabolic rate remained at near basal levels 

-2 -1 

and descended continuously from an initial 45 Kcal . m , h to a final 3 6 Kcal 
-2 -1 

. m . h . It will be seen from Table 15 that in the simulation some shivering 

started after about one hour in the cold, and it remained at one basal metabolic 

rate through much of the remaining hour. In the experiment, there was a con- 

-2 -1 

tinuing negative heat storage of about 37 Kcal . m , h throughout that same 

hour. In the experiment, it was calculated that the total net heat loss during 

-2 

the cold exposure amounted to 92 Kcal . m or about a 2.6°C drop in mean body 
temperature. We have also observed that the model response to this regimen of 
cold exposure is very close to that shown consistently by women subjects. 
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Figure 10 . Experiment in which subject spends 30 minutes at thermally 
neutral temperature; 120 minutes at ambient temperature of 
48° C; and 60 minutes recovery in 30 ° C environment. 





Figure 11. Experiment results 
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Performance of Model of Thermoregulation in Man During Exercise 


In order to test the performance of the model in predicting the results of ther- 
moregulatory stress resulting from internally produced heat during exercise, it was 
applied to conditions for which experimental measurements were obtained. Mini- 
mally dressed male subjects were seated on a bicycle ergometer, in an environ- 
ment of 30°C at 30% relative humidity, with an air velocity of 10 cm . sec \ 
Evaporative heat loss was determined by weight recording, and metabolic rate 
through recording of oxygen consumption. Figure 12 shows the experimental re- 
sults in solid lines, and the model predictions in broken lines. A more complete 
tabulation of the model responses is given in Table 16. This experiment constitutes 
a severe stress, and it requires well-trained individuals. Physically less fit sub- 
jects usually have to terminate the experiment before the end of the heavy exercise. 
During the final part of the heaviest exercise bout, the model estimates skin blood 
flow at over 3 liters . min \ and the percentage wetted area goes up to 66%, both 
these values are indices of thermal discomfort. 

A comparison of an identical exercise experiment carried out at an ambient 
temperature of 20°C and the model prediction for the same condition is shown in 
Figure 13, with more complete model predictions in Table 17. 

Discussion of Model Performance, Short Comings and the Path for Improvement 

Any model is by definition simpler than the system it attempts to represent. 
The structure and definition of a model is focused on certain characteristics or 
conditions which are deemed of primary importance. In the manned space program, 
the thermophysiological problems are mainly related to prevention of hyperthermia 
(Stolwijk (41) ). Both the experimental efforts and the emphasis of the model re- 
ported on here were directed at heat elimination. As can be seen in the previous 
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section, the model predicts with reasonable accuracy the dynamic thermoregulatory 
responses to dynamic loads of ambient temperature and internal heat production 
even during very heavy exercise. There are discrepancies of absolute value and 
sometimes for short periods even of direction of change. Such discrepancies 
are associated with two types of conditions: exposure to thermal environments 
in which cooling occurs, and during the onset of exercise. During prolonged 
cooling, such as shown in Figure 11, between 60 and 180 minutes the core cools 
more in the simulation than in the experiment, and consequently shivering is 
initiated. The shivering causes perfusion of cold muscle with warm blood, and 
a resultant further drop of core temperatures increases shivering even more. 

Similarly in Figure 13 at 3 5 minutes and again at 95 minutes, there is a marked 
drop of core temperature at the onset of exercise, also due to perfusion of 
relatively cold muscle with warmer blood. 

Since each of the segments in the controlled system consists of only 4 layers, 
there will be errors in the finite difference method of solution during the develop- 
ment of new gradients in the relatively thick muscle and core layers. During heat 
stress, such errors are small because the gradients are relatively small, and con- 
vective heat flow by the circulatory system is the major avenue of heat flow. During 
cooling this convective heat transport largely disappears, gradients appear and 
considerable errors can occur which tend to draw heat from the core in the model, 
before this happens in the experiment. Such errors tend to become smaller again as 
a steady state is approached. It is obvious that a simple remedy is available which 
consists of introducing additional layers in the controlled system. This method 
was used by Wissler (42) and it provides improved performance during cold exposure. 
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The erroneous, rapid and transient drop in internal temperature at the onset 
of exercise in the model, especially after prolonged rest in a cool environment, 
is probably the result of another simplification: immediately at the onset of 
exercise, the whole lumped mass of muscle in a given compartment is supplied 
with all the blood required for oxygen transport. Two factors are probably 
associated with the sudden drop in internal temperature. The present model 
does not make allowance for the development or repayment of an oxygen debt, and 
if such allowance were made it would tend to reduce the initial cooling of the blood 
in the muscle. A second factor might be that in a given muscle compartment not 
all muscles are likely to be active, so that those groups which are active warm 
up faster and cause less initial blood cooling than now occurs in the model. 

Since it is intended to combine the present model of thermoregulation with 
simultaneous models of the cardio-vascular and the respiratory system, it is to 
be expected that improvements in this particular idiosyncrasy of the present thermo- 
regulation model will result. 
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Table 16 


HlDtDlDtOUHOHHHfAJ-lDJtOjJOlJ-OrsN^H^tALOlAHOCICVjr^CMNCNOOtDiriJ-mCOMOl 



2<NCNfsl(NCNCNJHHHHNOJ-rACs|HHHOHOOMJtl£)HOOONinLPi H O O) J IA ID O H in CN H CO 
O- rH HHHHHHHCNCNrArAKUArvjHHHHKUALAtOiT^UDfAMCMCMH 


lAi-HOCDcnOTHi^r'-cMrALncNOi'^rAcoTArALArH.drAtDootDLDLnrHOCMLALA-d.d.dLDO.dtOLAP'. 



SrHrHrHOOiHLOLALALD-dOOCMCMrHrHOOCDLnLAO-dLOCDCMOCnOOP'-P^CNlLDIArs.OOO-dCDLO-dCM 

OHHHHHHHHHHCMCM(NfMCMCNICMMCMJ-LAOLOLOfArArA{NCMCNl^COCT)OOHCJ-infArAfA 

o HHH 


O 

O 


CONlfltDOOlCNHGNOOoOaiCOMDlAlDNOOinNlAJMMACMCniA^CTHLnj-inJCNOKlO 

lD(jDlDl01DL3J'J-KlfALnr'N(J)OOOOOOCOCOOOJfl s 'CnOHlAHOOOC^LA4-COH(N^Cnr^CNHHH 


j-^j-j-j-jr^NNf^Nr^j-j-^d-j-^OHHHc^csjiALAiniALAj'iDNrsoocoooLriuMnLnLALri 


HHHHHK\ClCnOiniA 00 im^NCMHOM|AK\OlAOaUDHcil 04 -aiNrAmcr. OOOOOUDHMA 
LLHHHHHHOOHH^J-t^lAKWOrONUnHJ-tDNcOMAlAJ^^tAJ-CnHN^jNaDlOmiA 

CC ............ 

COOOOOOOOOOOOOOOOOOOOrHrHrHrHrHOOOOOOrHCMCMIAtATArMOOOOO 


OHHNtAj-oooi^LOMHCsjrsm^i/untoocD^N^cocii^ooiorAj-oDinomoiooococoH 

4- (D N J-NOQ00N00HK>^d-J-H00lI3inrriCMt001OHCMCM00LACvjO00N 


l-LruniALrunmiDNNrNfN.SNSNioioioNoooocooocoooNNNNScocoaioiCTaioooooooosrN 

lArAmrAfnhAr^fA^r^f^rAr^mKunrAnfArAfArAf'UOrAfAfArAnfAfArAK'irAr^K^borArArArAn 


cn 

h 


CO 

h 


£o 


> 

UI 




b 


-d 

rs 

rH 

-d 

P- 

OO 

A* 

IA 

00 

OO 

O 

to 

CO 


to 

LA 

LA 

-d 

O 

r-H 

d 

TA 

OO 

o 

CO 

CD 

OO 

r*- 

to 

to 

CO 

TA 

rH 

O 

LA 

CD 

CM 

CD 

d 

CM 

rH 

O 

CD 

CD 

o 

O 

O 

O 

LO 

1^. 

00 

O 

CM 

CM 

rH 

r-H 

rH 

rH 

t-H 

t-H 

CM 

d 

LA 

to 

to 

p- 

CM 

rH 

r-H 

t-H 

rH 

»H 

LA 

CD 

rH 

CM 

CM 

CM 

LA 

CM 

CM 

CM 

CM 

CM 

CO 

(O 

p- 


p- 

K 

to 

to 

to 

P^ 


r- 


r-* 

P«* 


P^ 



A- 

P^ 

p- 

P^ 


P- 


P* 

P^ 

P- 

P^ 

P^ 

P^ 

00 

CO 

CO 

OO 


P- 

P- 

P- 

P^ 

P- 

TA 

TA 

TA 

TA 

TA 

|a 

TA 

TA 

TA 

TA 

IA 

TA 

TA 

TA 

TA 

TA 

TA 

TA 

TA 

TA 

TA 

TA 

TA 

TA 

TA 

IA 

TA 

IA 

IA 

TA 

rA 

TA 

TA 

IA 

IA 

IA 

IA 

IA 

IA 

TA 

IA 

A 

00 


LA 

00 

o 

o 

CD 

CD 

LO 

d 

CM 

to 

00 

P^ 


to 

to 

LA 

CD 

d 

d 

O 

-d 

to 

CD 

CD 

OO 

P^ 

P^ 

to 

CD 

rH 

CM 

00 

CM 

LA 

OO 

d- 

A 

CM 

rH 

O 

P^ 

°o 

OO 

OO 

CD 

p 

-3- 

to 

OO 

O 

r-H 

«— 1 

CD 

CD 

CD 

CD 

CD 

CD 

CM 

d 

LA 

to 

to 

to 

CD 

CD 

CD 

CD 

CD 

CD 

0O 

rH 

CM 

CM 

IA 

TA 

o 

O 

O 

o 

O 

O 

to 

LO 

to 

to 

to 

Lo 

to 

to 

to 

A- 

Ps 

b 

LO 

to 

to 

to 

to 

to 

p* 

fs. 

P^ 

p^ 


p^ 

to 

to 

to 

to 

to 

to 

p" 

CO 

00 

00 

OO 

OO 

Ps 

r- 


p- 

P^ 

P- 

IA 

A 

TA 

TA 

TA 

hA 

TA 

IA 

IA 

TA 

IA 

IA 

TA 

TA 

IA 


TA 

TA 

TA 

TA 

IA 

IA 

IA 

TA 

TA 

TA 

IA 

IA 

IA 

IA 

TA 

TA 

TA 

TA 

IA 

TA 

IA 

IA 

IA 

rA 

TA 

IA 

OO 


tA 

OO 

»-H 

1m 


u 


P* 

00 

Ld 

CM 

rH 

o 

CD 

CD 

CD 

d 

LO 

O 

CD 

-d 

P* 

to 

J- 

CM 

rH 

rH 

O 

o 

CO 

Pv 

P^ 

IA 

LD 

CM 

LA 

00 

to 

LA 

d 

CD 

i° 

O 

O 

rH 

p 


CO 

CD 

r-H 

CM 

P 

CM 

CM 

CM 

r-H 

rH 

M 

IA 

LA 

P>- 

PS 

00 

00 

TA 

CM 

CM 

CM 

CNI 

CM 

00 

|H 

IA 

d 

LA 

LA 

LO 

IA 

CM 

CM 

CM 

CM 

to 

K 

p- 

p^ 

Ps 

K 

LO 

to 

LO 

p- 

p^ 

P 


P^ 

rs. 

pT 

P^ 

P 


P^ 

P^ 


p^ 


P^ 

P^ 

P^ 

P^ 

P^ 

r>- 


OO 

CO 

oo 

OO 

CO 

P- 

P>- 

Ps 

P- 

P- 

P^ 

IA 

ta 

IA 

TA 

TA 

h 

TA 

h 

TA 

TA 

IA 

TA 

j 

IA 

TA 

IA 

n 

TA 

fA 

IA 

fA 

TA 

TA 

TA 

TA 

IA 

IA 

IA 

IA 

IA 

TA 

IA 

TA 

IA 

TA 

IA 

TA 

TA 

IA 

IA 

TA 

IA 

IA 

P^ 

b 

Ps 

jp* 


b 

r-H 

p 

o 

LA 

r-H 

u 

CD 

id 

LA 

TA 

CM 

b 

CD 

f 

o 

O 

to 

CD 

rH 

CD 

to 

LO 


rH 

p^ 

r~- 

O 

CM 

CD 

LD 

rH 

CM 

O 

rH 

OO 

d 

O 

CD 

p 

CD 

CD 

CD 

P 

00 

p 

o 

H 

CM 

Km 

j . 

-d 


LA 

LA 

LA 

P 

TA 

p 

p-. 

p- 

p- 

CO 

p* 

OO 

0O 

CO 

oo 

h- 

CO 

rH 

CM 

CM 

TA 

d 

O 

O 

O 

CD 

CD 

CD 

rA 

(A 

TA 

TA 

TA 

k 

IA 

b 

d 

d" 

.d 

Ld 

d" 

-d 

-d 

hd 

d 

id 

d 

d 

-d 

-d 

d 

d 

-d 

d 

-d 

J- 

-d 

d- 

-d 

LA 

LA 

LA 

LA 

LA 

LA 

LA 

LA 

d- 

d 

-d 

TA 

fA 

TA 

TA 

TA 

n 

TA 

IA 

TA 

TA 

TA 

TA 

1 

TA 

TA 

TA 


IA 

TA 

TA 

TA 

TA 

TA 

IA 

TA 

IA 

IA 

IA 

IA 

TA 

TA 

TA 

IA 

IA 

IA 

IA 

TA 

TA 

TA 

TA 

ta 

TA 

TA 

to 

CO 

CD 

O 

t — 1 

i 

pj 

iH 

^d 

o 

CD 

IA 

^A 

LA 

rH 

to 

CM 

OO 

d 

CM 


-d 

LA 

TA 

P^ 

LA 

TA 

to 

CD 

CM 

LO 

to 

CO 

P^ 

Ps 

d- 

O 

CM 

d 

rH 

rH 

t — 1 

CM 


h 

P^ 

OO 

OO 

po 

O 

CM 

d 

LA 

rs. 

CO 

r- 

r* 

to 

to 

LA 

in 

CD 

rH 

IA 

d 

LA 

LA 

IA 

CM 

t—l 

O 

O 

CD 

LA 

CO 

O 

rH 

CM 

TA 

CD 

P- 

LD 

LA 

d 

TA 

LA 

IA 

LA 

LA 

LA 

LA 

LO 

to 

to 

to 

to 

P 

to 

to 

to 

to 

to 

to 

to 

P-I 

P^ 

P- 

P^ 

P^ 


P^ 

P^ 

P^ 

P^ 

to 

P~«* 

p**- 

OO 

oo 

OO 

OO 

P** 


P- 

P^ 

P^ 

P- 

TA 

m 

TA 

(A 

TA 

TA 

TA 

IA 

TA 

TA 

TA 

fn 

IA 

TA 

IA 

TA 

IA 

TA 

TA 

IA 

IA 

TA 

TA 

TA 

IA 

IA 

TA 

IA 

IA 

TA 

TA 

IA 

IA 

TA 

IA 

A 

TA 

TA 

TA 

TA 

A 

IA 


U5l01DUMDlO»^lArOCOtnpl^KUAtnNOOCMJ-OUASrvCTlNh*Nl^(Nd-oONHCDXhOOCSncn 


IS n n n n j-j-j-JtNcn^inj-rAMMOooiDOHNrscoj-ojOcncMLDLri j- m rs 01 oo oo rA o p* 

HHHHrnj-NCNCMCvjCSHNiriOCMJ'lALAinKUnnrACvjLnCMd) OO CD CD o m J- d d !A 

rH r-H r-H i — I I — I tH CM CM CM CM CM «H 




OOlCDOlOIONMA^^^aiaiOlOlOlOwrinKUAIA^aiCJlOlClOlCnOOOOOOOlOIOHJlOIOl 
-diAiATAiATAcMCMCMCMCMCMrArArAiAiATACMCMCMcMCMtMTAiAiATAiATAP^r^p'.p-p^p-iATATATATATA 
HHH HHH CM M CM (N CM N IA IA IA |A TA TA 


HCOLnj-^HlD^lCMtACTiKiMmAJ-lACONj-OCPOin^OCONlAOfAcOOJjtHINJ-MHO^ 


LA 

d 

d 

d- 

-d 

-d 

d 

IA 

o 

P- 

TA 

rH 

CD 

P- 

LO 

LA 

d 

TA 

O LD 

CD 

TA 

TA 

to 

rH 

TA 

00 

LO 

d 

IA 

CO 

iH 

rH 

rH 

CM 

P- 

LA 

TA 

TA 

CO 

LA 

CM 







PS 

Ps 

p- 

to 

d 

TA 

rH 

rH 

rH 

rH 

rH 

rH 

CM 

r- 

-d 

IA 

CM 

rH 

LA 

IA 

CM 

CM 

CM 

CM 

to 

CD 

LA 

TA 

CM 

rH 

O 

LA 

d 

IA 

TA 

TA 













1 

1 

' 

J 

1 

1 

t — 1 






1 

1 

1 

1 

1 

' 

rH 






i — 1 

1 

1 

1 

1 

1 

LA 

O 

LA 

O 

LA 

O 

LA 

O 

LA 

o 

LA 

a 

LA 

O 

LA 

O 

LA 

o 

LA 

o 

LA 

O 

LA 

O 

LA 

O 

LA 

O 

LA 

o 

LA 

O 

LA 

O 

LA 

O 

LA 

O 

LA 

O 

LA 

o 


rH 

rH 

CM 

CM 

TA 

TA 

-d 

d 

LA 

LA 

LD 

to 

P^ 

P^ 

CO 

OO 

CD 

CD 

o 

O 

rH 

rH 

CM 

CM 

IA 

IA 

d 

d 

LA 

LA 

to 

to 

P- 

P- 

OO 

CO 

CD 

CD 

o 

O 

rH 




















rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

1 — 1 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

CM 

CM 

CM 


69 







Table 17 


hNMCMCMMMOOOOOOCMMMCNMCNOOOOj-CJlNtDtNNNMO^molOOlArNHOOHOfN 

• • 

^CMfSICMCM(NCMHHHHHH(NCMCM(MCNCNlHHHHCOr^OCVJMCMCMCNHCDLnrAoOOHj-OOLOJ-(N 
0- rH «H rH K\ J- ^ lA ^ H 


HJ > fAH^NNOOJ , CNOCTllANlDJ , NOMOCT)H 4 -CNKOK\HlOOtAOHHl/UDCOr><.inOJ , NM 



ZHrsLnj-fACMQinirn^iAJ-HHHHHHiAiA^inaiMLaj'^rru^Mtomj-HoiNooHai^txjin 

CCNHHHHHHHHHHHHHHHHHHHHHHCSHHHHHHHCvjfAJ-J-lACNJCNHHHH 


ClHlJLANHOOOfs/l^incnCNLOJ'lAOOlAJ'CTOOrAinrAHtDJ'CNHOJ-LnOCDOOOJ-pNCNCnNa 

(£)NNoOCnHLOinMOCOtD01CCOOOOCOaHNlOI^4lAlDNLDCDlfltDlDOl^N01NlDOlNNtDlDlC 



U^44J^^00C00003NN^JtJJ44HOOOOO4 4-14J^ininimAl0l0J4JJ4^ 



CO 

LA 

.3 

st 

st 

A A 

A 

A A 

A 

A 

A 

A 

A A 

A 

A A 

A 

A CO O 

00 CO rH 

07 p 

CO 

A 

CO 

07 

St CM 

CM 

st OO 

1 tH 

CO -zf 

CJ 

tH 

LL 

O 

o 

O 

o 

O 

O O 

O 

O O 

O 

O 

O 

O 

o o 

O 

o o 

O 

O O CM 

CM H rH 

o o 

o 

O 

O 

A 

o 

A 

CD A 

t CM 

rH rH 

rH 

rH 

00 

o 

o 

O 

o 

O 

o o 

O 

O O 

O 

O 

o 

O 

o o 

o 

o o 

O 

o o O 

0. 

0. 

0. 

o o 

o 

o 

o 

O 

O rH 

rH 

rH O 

o 

o o 

O 

O 


CO 

LA 

rH 

CD 

p^ 

CO -3 

-3 

CD CM 

CO 

CO 

to 

07 

st cm 

A 

A ID 

07 

A st A 

p a o 

p o 

A 

A 

o 

CD 

07 CM 

A 

o o 

00 

O 07 

CM 

CD 


CM 

rH 

O 

oo 

p* 

CO oo 

A 

A P- 

CO 

07 

st 

O 

oo to 

St 

CM CO 

A 

nom 

A 00 st 

O OO 

A 

A 

CM 

P» 

O CM 

CM 

CM P^ 

CM 

07 A 

A 

O 

1— 

A 

LA 

A 

st 

st 

st A 

CD 

CO CO 

CO 

CO 

to 

to 

A A 

A 

A CD 

p 

P CO CO 

00 N S 

p to 

CO 

CD 

00 

00 

CD 07 

CD 

07 00 

oo 

P^ P^ 

P^ 

P^ 


A 

A 

A 

A 

A 

A A 

A 

A A 

A 

A 

A 

A 

A A 

A 

A A 

A 

AAA 

AAA 

A A 

A 

A 

A 

A 

A A 

A 

A A 

A 

A A 

A 

A 


a 

St 

CM 

OO 

rH 

-3 p- 

rH 

cd 

rH 

rH 

rH 

A 

CO A 

A 

O rH 

07 

A St St 

CO st CO 

A rH 

rH 

O 

07 

07 

st O 

A 

rH CD 

CO 

St A 

CM 

CM 


CD 

CD 

CD 

00 

00 

cm 

A 

St A 

P 

CO 

07 

CD 

CD CD 

07 

CD A 

P 

H st CO 

p a st 

St st 

St 

st 

CM 

P^ 

tH A 


A A 

st 

St St 

St 

St 

1- 

CO 

CO 

CO 

CO 

CO 

CO CO 

CO 

CO CO 

CO 

to 

CO 

to 

CD tO 

CO 

CO CO 

CO 

p p p 

p p p 

p p 

p 

p 

p 


oo oo 

00 

OO P-* 

p^ 


p^ 

P^ 


a 

A 

A 

A 

A 

A A 

A 

A A 

A 

A 

A 

A 

A A 

A 

A A 

A 

AAA 

AAA 

A A 

A 

A 

A 

A 

A A 

A 

A A 

A 

A A 

A 

A 


CO 

CO 

CM 

CO 

CD 

rH OO 

IP 

A -3 

A 

-3 

A 

P 

00 to 

St 

O CO 

A 

00 st CO 

P OO A 

A CM 

CM 

CM 

A 

rH 

A A 

CM 

A tO 

St 

st st 

tn 

CM 


p- 


p^ 

CO 

A 

A O 

CM 

-3 A 

CO 

ip 

p 

p 

P p 

P 

P A 

07 

CM A CO 

P CM CM 

CM CM 

CM 

CM 

CO 

iH 

A st 

st 

A CM 

CM 

CM CM 


CM 

1- 

CO 

CO 

CO 

CD 

CO 

CO CO 

CO 

CO CO 

CO 

to 

CD 

CO 

to to 

CO 

CO CO 

CO 

p p p 

p p p 

p p 

P 

P 

p 

OO 

CO oo 

<P0 

00 Ps. 

P^ 

P- P-. 

k‘ 

P» 


A 

(A 

A 

A 

A 

A A 

A 

A A 

A 

A 

A 

A 

A A 

A 

A A 

A 

AAA 

AAA 

A A 

A 

A 

A 

A 

A A 

A 

A A 

A 

A A 

A 

A 


p* 

CO 

CM 


O 

CM CM 

CO 

rH A 

St 

A 

CO 

O 

rH rH 

07 

CO O 

CD 

rH o O 

CM H CM 

07 OO 

P 

to 

A 

CO 

A CM 

P 

p** to 

CD 

'st CM 

rn 

O 


CD 

CD 

CD 

00 

OO 

p. A 

A 

A CO 


00 

00 

07 

CD 07 

00 

OO to 

00 

Nins 

00 A st 

A A 

A 

A 

ft 

CD 

A A 

CO 

A tO 

-zf 

r* st 

ft 

St 

i— 

CO 

CO 

CO 

CO 

CO 

CO CO 

CO 

CO CD 

to 

CO 

to 

to 

to to 

to 

CO tO 

to 

p p P 

P P P 

p p 

P 

P 

K 

P^ 

<po 00 

to 

OO P- 

P^ 

P- p*. 

S 

p. 


a 

A 

A 

A 

A 

A A 

A 

A A 

A 

A 

A 

A 

A A 

A 

A A 

A 

AAA 

AAA 

A A 

A 

A 

A 

A 

p A 

A 

A A 

A 

A A 

h 

A 


oo 

st 

CO 

^t 

CO 

A O 

CD 

CD 07 

rH 

rH 

rH 

p 

A tO 

OO 

07 rH 

OO 

St p St 

A O -3 

OO st 

CO 

rH 

P^ 

CM 

CD 


O CM 

A 

00 A 

b 

P^ 


st 

oo 

A 

o 

IP 

A CD 

rp 

P* IP 

CO 

CO 

A 

A 

A CM 

rH 

O OO 

CD 

rH A CO 

P A -3 

A CM 

P 

07 

CO 

O 

A C- 

r 

OO -zf 

CM 

O 07 

p 

CO 

l- 

CM 

rH 

rH 

rH 

o 

O CD 

CD 

07 07 

CD 

07 

O 

o 

OO 

o 

O 07 

07 

ooil 

O rH rH 

iH rH 

tH 

O 

o 

iH 

H rH 


CM 

CM 

[ • • 

fcM rH 

M 

rH 


r A 

A 

A 

A 

A 

A CM 

CM 

CM CM 

CM 

CM 

A 

A 

A A 

A 

A CM 

CM 

AAA 

AAA 

A A 

A 

A 

A 

A 

A A 

A 

A A 

A 

A A 

^7 

A 


o 

fp 

st 

A 

CO 

CD CD 

st 

00 CO 

A 

CM 

A 

P 

07 rH 

A 

A OO 

07 

CO rH A 

N07 tn 

CM O 

CO 

CO 

±t 

A 

CO OO 

CM 

P- P^ 

A 

CD rH 

P^ 

A 

rr> 

CO 

st 

A 

CM 

rH 

O rH 

A 

-3 A 

CO 

p 

CO 

A 

st St 

A 

CM A 

OO 

rH -3 CO 

P A st 

A CM 

o 

07 

A 

O 

A A 

N 

P^ A 

rH 

CD 00 

b 

A 

1— 

cn 

A 

A 

A 

A 

A A 

A 

A A 

A 

A 

A 

A 

A A 

A 

A A 

A 

CO CO to 

CO CO CO 

CO CO 

to 

A 

CO 

P^ 

r>! p^ 


P-T p! 


CO CD 

to 

CO 


a 

A 

A 

A 

A 

A A 

A 

A A 

A 

A 

A 

A 

A A 

A 

A A 

A 

AAA 

AAA 

A A 

A 

A 

A 

A 

A A 

A 

A A 

A 

A A 

A 

A 


st 

A 

CM 

rH 

O 

O -zi- 

-3 

-3 3- 

-3 

st 

O 

O 

o o 

07 

07 CM 

A 

A st O 

3 

5 

8 

CM CM 

CM 

i — 4 

A 

P^ 

rH LO 

rH 

O CM 

O0 

CM A 

rH 

A 

> 

c- 


P- 

p^ 

|p 

p'^ -3 


st -3 

st 

-3 

p 

p 

p p 

CO 

CO A 

A 

A A OO 

CO 07 N 

P P 

p 

P 

CO 

OO 

• « 

A OO 

• 

rH 

• • 

A CD 

• 

A 

• • 

OO A 

O 

• 

LU 






rH 

rH 

rH rH 

rH 

rH 





CM 

CM 

CM CM _3 

tO rH 




A 

07 

A CO 

rH 

CM A 

CM 

rH rH 

rH 


























rH rH 

CM 

CM 






CD 

oo 

p. 

p* 

CO 

CM rH 

A 

rH O 

A 

CM 

rH 

A 

A P 

OO 

A CM 

A 

AAA 

A st st 

st St 

st 

st 

CM 

CM 

CM CM 

CM 

CM St 

st 

-Zf -zf 

*3- 

-3 

2 : 

i— 1 

CM 

-3 

r^. 

rH 

CO P" 

CO 

CO O 

A 

OO 

O 

oo 

p p 

oo 

O rH 

A 

AAA 

A 07 07 

07 07 

CD 

07 

O 

O 

O O 

O 

O CD 

07 

07 07 

' • 

CD 

07 


-3 

J- 

st 

st 

A 

A CO 

CO 

A A 

-3 

A 

A 

st 

St St 

St 

A A 

A 

CM CM CM 

CM A A 

A A 

A 

A 

P-* 


P- 

P-* 

P-* A 

A 

A A 

A 

A 







rH 

rH 

rH rH 

rH 

rH 





CM 

CM 

CM CM CM 

CM 




A 

A 

A A 

A 

A 





CD 

A 

CD 

A 

|p 

A CM 

CD 

CM 07 

rH 

00 

CD 

A 

oo rH 

St 

iH CM 

to 

8 

8 

9 

3 
0 

4 

A A 

A 

rH 

rH 

o 

P- O 

OO 

CD OO 

O 

A CO 

CM 

A 

CO 

rH 

CO 

o 

A 

CD 

A CM 

CM 

St A 

07 

A 

p 

O 

O O 

oo 

CO oo 

CO 

CD St CO 

n oo rs 

CO A 

St 

A 

A 

07 

0O rH 

A 

• • 

CO CM 

• 

CM 

• • 

CM tO 

CM 

00 


a 

St 

.3 

A 

CM 

CM A 

A 

-3 A 

CM 

CM 

CM 

A 

A A 

CM 

CM CM 

o 

07 07 tO 

st A st 

St St 

st 

-3 

CM 

A 

07 CO 

A 

rH O 

P. 

CO A 

A 

-3 ' 


i 

1 

1 

1 

1 

1 





' 

1 

1 1 

1 

1 rH 

rH 


1 1 

1 1 

1 

t 

CM 

rH 



rH 

1 

1 

1 1 

* 

1 : 

UJ 

LA 

o 

A 

O 

A 

O A 

O 

A O 

A 

O 

A 

O 

A O 

A 

O A 

O 

A O A 

O A O 

A O 

A 

O 

A 

O 

A O 

A 

O A 

O 

A O 

A 

o 

— 


rH 

rH 

CM 

CM 

A A 

st 

st A 

A 

to 

CO 

p 

p oo 

OO 

07 07 

O 

O rH rH 

CM CM A 

A st 

St 

A 

A 

CO 

to P-* 

P» 

CO oo 

07 

07 O 

O 

rH 

h- 
















rH 

rH rH rH 

H H H 

rH rH 

rH 

rH 

rH 

rH 

rH rH 

rH 

rH rH 

rH 

rH CM 

CM 

CM 


71 


REFERENCES 


Stolwijk, J. A. J. and J. D. Hardy, Temperature regulation in man - 
a theoretical study. Pflugers Arch. 291 , 129-162 (1966). 

Yamamoto, W. S. and W. F. Raub, Models of the regulation of 
internal respiration in mammals. Problems and promises. 
Computers and Biomet. Res. _1, 65-104 (1967). 

Meeh, K. , Oberflachenmessungen des menschlichen Korpers, J. Biol., 
15, 428-458 (1879). 

DuBois, D. and E. F. DuBois, Clinical Calorimetry: Fifth paper. The 
measurement of surface area of man. Arch. Intern. Med., 

15 , 868-888 ( 1915 ). 

Bradfield, H. S., The determination of surface area of women and its 
use in expressing basal metabolic rate. Am. J. Physiol., 

82, 570-576 (1927). 

Wilmer, H. A., Changes in structural components of human body from 
six lunar months to maturity. Proc. Soc. Exp. Biol. Med., 

43, 545-547 (1940). 

Scammon, R. E., Developmental Anatomy In Schaeffer, J. P. (Editor) 
Morris' Human Anatomy, 11th Ed., New York, McGraw Hill 
(1953). 

DuBois, D. and E. F. DuBois, Clinical Calorimetry: Tenth paper, 

A formula to estimate the approximate surface area if height 
and weight be known. Arch. Intern. Med., _17, 863-871 (1915). 



REFERENCES 


9. Nishi, Y. and A. P. Gagge, Mean skin temperature for man weighted 

by regional heat transfer coefficients. Fed. Proc. 19, 852 
(1970). 

10. Stolwijk, J. A. J. and J. D. Hardy, Partitional calorimetric studies of 

responses of man to thermal transients. J. Appl. Physiol., 

21 , 967-977 (1966). 

11. Hardy, J. D. and J. A. J. Stolwijk, Partitional calorimetric studies of 

man during exposures to thermal transients. J. Appl. Physiol., 

21 , 1799-1806 (1966). 

12. Brebner, D. F., D. McK. Kerslake and J. L. Waddell, The relation 

between the coefficients for heat exchange by convection 
and evaporation by man. J. Physiol., 141 . 164-168 (1958). 

13. Gagge, A. P., A new physiological variable associated with sensible 

and insensible perspiration. Am. J. Physiol., 120 , 277-287 
(1937). 

14. Gagge, A. P., J. A. J. Stolwijk and B. Saltin, Comfort and thermal 

sensations and associated physiological responses during 
exercise at various ambient temperatures. Environmental Research, 

2, 209-229 (1969). 

15. Aschoff, Y. and R. Wever, Kern und Schale im Warmehaushalt des Meuschen. 

Naturwinenschaften, .45, 477-485 (1958). 

16. Scheinburg, P., L. I. Blackburn, M. Rich and M. Saslaw, Effects of 

vigorous physical exercise on cerebral circulation and 
metabolism. Am. J. Medicine, JJ5, 549-554 (1954). 


73 



REFERENCES 


17. Behnke, A. R. and T. L. Willmon , Cutaneous diffusion of helium in 

relation to peripheral blood flow and absorption of 
atmospheric nitrogen through the skin. Am. J. Physiol., 

131 . 627-632 (1941). 

18. Hertzman, A. B. and W. C. Randall, Regional differences in the basal 

and maximal rates of blood flow in the skin. J. Appl. Physiol., 

I, 234-241 (1948) 

19. Cabanac, M., J. A. J. Stolwijk and J. D. Hardy, Effect of temperature 

and pyrogens on single-unit activity in the rabbit's brain 
stem. J. Appl. Physiol., 24 . 645-652 (1968). 

20. Hellstr^m, B. and H. T. Hammel, Some characteristics of temperature 

regulation in unanesthetized dog. Am. J. Physiol., 

213 . 547-556 (1967) 

21. Simon, E., W. Rautenberg, R. Thauer and M. Iriki, Die Auslosung von 

Kaltezittern durch lokale Kuhlung im Wirbelkanal. Pflugers Archiv., 
281 , 309-331 (1964). 

li II 11 

22. Wunnenberg, W. and K. Brack, Untersuchimgen uber die Funktionsweise 

thermoreceptiver Strukturen im Cervicalmark des Meerschweinchens . 
Pflugers Archiv. , 294 , 83 (1967). 

23. Hensel, H., Electrophysiology of cutaneous thermoreceptors. In. 

Skin senses, D. R. Kenshalo, Editor, Thomas, Springfield, 111., 
(1968). 

24. Bullard, R. W. , M. R. Banerjee and B. A. MacIntyre, The role of the skin 

in negative feedback regulation of eccrine sweating. Int. J. 
Biometeor., H.i 93-104 (1967). 


74 



REFERENCES 


25. Rawson, R. O. and K. P. Quick, Evidence of deep-body thermoreceptor 

response to intra-abdominal heating of the ewe. J. Appl. 
Physiol., 28, 813-820 (1970). 

26. Wood, J. D., Electrical activity from single neurons in Auerbach's 

plexus. Am. J. Physiol., 219, 159-169 (1970). 

27. Robinson, S., F. R. Meyer, J. L. Newton, C. H. TS'ao and L. O. 

Holgersen, Relations between sweating, cutaneous blood 
flow and body temperature in work. J. Appl. Physiol., 20 , 
575-582 (1965). 

28. Bader, M. E. and M. B. Macht, Indirect peripheral vasodilatation 

produced by warming of various body areas. J. Appl. Physiol., 
I, 215-226 (1948). 

29. Randall, W. C., Quantitation and regional distribution of sweat glands 

in man. J. Clin. Invest., 25., 761-767 (1946). 

30. Stolwijk, J. A. J., B. Saltin and A. P. Gagge, Physiological factors 

associated with sweating during exercise. J. Aerospace Med., 
39, 1101-1105 (1968). 

31. Bishop, J. M., K. W. Donald, S. H. Taylor and P. N. Wormald, 

Changes in arterial hepatic venous oxygen content difference 
during and after supine leg exercise. J. Physiol., 137 . 309-317 
(1957). 

32. Bucht, H., J. Ek, H. Eliasch, A. Holmgren, B. Josephson and L. Werko, 

The effect of exercise in the recumbent position on the renal 
circulation and sodium excretion in normal individuals. Acta 
Physiol. Scand., 28, 95-100 (1953). 


75 



REFERENCES 


133 

33. Grimby, G. E., E. Haggendahl and B. Saltin, Local Xenon 

clearance from the quadriceps muscle during exercise 
in man. J. Appl. Physiol., 22^ 305-310 (1967). 

34. Saltin, B., G. Blomquist, J. H. Mitchell, R. Johnson, Jr., 

K. Wildenthal and C. D. Chapman, Response to submaximal 
and maximal exercise after bed rest and training. 

Circulation, 3J3, Suppl. 7, 1-78 (1968). 

35. Saltin, B., A. P. Gagge and J. A. J. Stolwijk, Muscle temperature 

during submaximal exercise in man. J. Appl. Physiol., 25_ , 

679-688 (1968). 

36. Webb-Peploe, M. M. and J. T. Shepherd. Peripheral mechanism 

involved in response of dogs' cutaneous veins to local tem- 
perature change. Circulation Res ., 23., 701-708 (1968). 

37. Wurster, R. D. and R. D. McCook, Influence of rate of change of 

skin temperature on sweating. J. Appl. Physiol. 2J., 237-240 (1969). 

38. Nadel, E. R., R. W. Bullard and J. A. J. Stolwijk, Sensitivity to rate 

of change in skin temperature in the regulation of sweating, 
m.s. to be submitted. (1970). 

39. Benzinger, T. H., C. Kitzinger and A. W. Pratt, The Human Thermostat 

in Temperature - its measurement and control in science and industry, 
vol. 3, J. D. Hardy, Editor, p 637-665, Reinhold, New York (1963). 


76 


REFERENCES 


40. Nadel, E. R., The role of peripheral thermoreceptors in the integrated 

thermoregulatory response in resting and exercising man. 
Dissertation, Univ. Calif., Santa Barbara (1969). 

41. Stolwijk, J. A. J., Thermal loads in lunar ambulation. Aerospace 

Medicine (in press) (1970). 

42. Wissler, E. H. f A mathematical model of the human thermal system. 

Bull. Math. Biophysics, 2_6, 147-166 (1964). 


NASA-Langley, 1971 4 CR— 1855 


77 



